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ABSTRACT 
This study developed and compared four econometric models 
for accuracy and ease of use in predicting milk deliveries 
in Alberta for the time period November 1978 to June 1979. 
Two of the models were designated as Complex models and were 
developed by initially segregating dairy herds into five 
unique producer categories and then using ordinary least 
squares techniques to estimate the parameters for each 
month’s predictive equations. One of these models was 
designated as a Complex Stage-of-Lactation model and was 
developed from data encompassed in the Cow Master file of 
the Dairy Herd Improvement (DHI) program. The other model 
was referred to as a Complex Cows-in-Milk model and was 
developed using data encompassed in the Herd Master file of 
the same program. The other two models were designated 
Simplified Stage-of-Lactation and Simplified Cows-in-Milk 
and were developed from a smaller sample base and 
irrespective of producer category. 

Monthly delivery schedules were compared over the 12 
month period July 1978 to June 1979 to allow herd 
comparisons between DHI and non DHI herds for each category. 
Monthly delivery and production totals over the same time 
period were compared to assess the proportion of milk 
produced that was actually delivered to dairy plants. 

Total milk delivered in Alberta for a month under 
investigation was estimated to be the product of (a) 


predicted production per category per month for DHI herds, 
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(b) the number of herds in each category, (c) the monthly 
production coefficient, and (d) the monthly herd delivery 
coefficient. 

The Complex models were found to predict deliveries to 
within 2% of actual milk delivered and the Simplified models 
had an average 5.5% predicted error when the two model 
designations were compared over the same time period used in 
their development. Due to present data limitations the 
Complex Cows-in-Milk model was chosen as most appropriate 
for further research and completion of the model for year 
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I. INTRODUCTION 


A. The Problem 

In recent years Alberta dairy policy makers have faced 
the problem of increasing Alberta's total milk production 
and increasing Alberta’s share of industrial milk produced 
at the national level. In 1978 this national share amounted 
to approximately 7%' but a buoyant provincial economy has 
made the goal of increased production and increased market 
share a major concern for the 1980s. 2 

Alberta dairy policy, whether from an explicit or 
implied point of view, has always been attempting to 
increase provincial milk production. Over the years these 
attempts have been made in the face of declining yearly 
production and increasing average production per cow.? Today 
cow numbers in Alberta number 142,000 (excluding dairy 
heifers!) /va drop of 44.trom that of 19/78 2Figure: 1. 1 
visually shows this continuing drop in cow numbers. With 
Alberta's population increasing at approximately 3%4 and per 
capita fluid milk consumption increasing by over 5%5 the 
search for a stimulus to increase dairy production to cover 


provincial consumption is felt to be strongly needed. 


' Dairy Farmers of Canada 1979 Dairy Facts and Figures at a 
Glance Ottawa May 1979 Table 23 

2 Alberta Agriculture Alberta Dairy Situation and Qut look 
Market Analysis Analysis Branch April 1979 pt 

3 Ibid, Table 8 and Table 5 

4 Agriculture Canada Handbook of Food Expenditures, Prices 
and Consumption Ottawa April 1979 Table 2 pp 5 6 

5 Alberta Agriculture Alberta Dairy Situation and Qutlook 
Market Analysis Branch February 1979 Table 4 p3 
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B. Characteristics of the Alberta Dairy Industry 

Although the dairy industry in Alberta does not play a 
major role in the structure of the Alberta agricultural 
scene it is one industry that offers highly stable 
production of an important commodity. In 1976 Alberta had 
1745 commercial farmers who produced 581.3 thousand tonnes 
of milk.® By 1978 production had fallen to 577.3 million 
tonnes after reaching 584.7 tonnes in 1977.7 Thus milk 
production fluctuated little over the 3 year period. On a 
farm receipt basis, Alberta’s dairymen received 8% of the 
cash receipts of Alberta producers for livestock and 
livestock products and 5% of the total farm cash receipts of 
Alberta agricultural producers. 8 

The increase in the dairy heifer population to 31,000 
(a 6.9% increase from that of 1978) should have been an 
encouraging sign for dairy policy planners. However, in 
spite of this increase, the milk cow-heifer ratio rose to 
4.6 versus the ratios of 4.3 and 4.1 for 1977 and 1978 
respectively. This increase in the ratio suggests a 
depletion in the number of replacement cows and could be 
signalling future declines in both dairy cow numbers and 
milk production. As well this higher ratio could be due to 
rising feeder cattle prices which makes more stringent herd 


6 Dairy Farmers of Canada 1979 Dairy Facts and Figures at a 
Glance Ottawa May 1979 Table 6 

cared 

8 Statistics Canada Farm Cash Receipts Ottawa February 1979 
Vol 39 No 12 
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culling and the raising of dairy heifers for beef purposes 
financially feasible. To counteract these possible 
detrimental circumstances a number of recent policy changes 
have been made.’ Unfortunately it is difficult to assess the 
dairy industry’s response to the changes because there is no 
reliable system of measurement to assess dairymen’s plans 
for expanding herd sizes, increasing average herd production 
or both. 

This study will attempt to overcome this problem. Data 
from the Dairy Herd Improvement Program will be used as a 
base for the development of four econometric milk predicting 
models. Multiple regression or ordinary least squares 
analysis will then be used to estimate the coefficients of 
each model. Since not all herds in Alberta are enrolled on 
the program, a comparison of milk deliveries to respective 
plants will be undertaken to compare DHI and non-DHI 
producers for a common type of shipment. This data is 
compiled monthly by the Alberta Dairy Control Board. The 
proportion of milk produced that was actually delivered will 
be determined since the supply of milk in Alberta is 
measured by the total deliveries made to dairy plants By 
accurately predicting milk supply, policy planners should be 
better able to anticipate the industries response to present 
dairy policy and make the necessary changes to encourage 


production. 


9 These policy changes will be discussed more fully later 
on. 
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The price of milk and the prices of other complementary 
and competing products will not be used in the development 
of the models because the dairy industry is highly 
regulated. Any changes that occurred in these areas in the 
time period under study will be viewed as being "stochastic" 
in nature. These models will be an attempt to measure 
industry response after it has occurred as opposed to a 
larger and more complex supply response model which usually 
attempts to measure industry response in relation to a 


change in one variable, usually the price of milk. 


C. Limitations of the Study 

This study will be done using data over a short time 
period. Since DHI data is not stored on a regular basis, the 
data unique to development of two of the four models could 
only be obtained on an “as it happens" monthly basis. Hence, 
DHI data for only the months of November 1978 to June 1979 
inclusive were used in this study. Selection of the most 
appropriate model for continued research will be made from a 


comparison of the models over this time period. 


D. Organization of the Thesis 

A discussion of dairy policy at both the Federal and 
Provincial level follows in the next chapter. It also 
contains an overview of the changes that have been made to 
current policy and the economic theory that serves as a 


background to these changes. Chapter 3 contains a review of 
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previous research work done in the area of predicting milk 
production. An outline of the research methodology used in 
the development of the four models is the topic of 
discussion in Chapter 4 and and Chapter 5 contains 
discussion on the results, accuracy and suitability of the 
four models involved. The final chapter includes a brief 
summary of the study, the conclusions drawn and a discussion 


on the directions for future research. 
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II. AN OVERVIEW OF GOVERNMENT DAIRY POLICY 

Canadian dairy policy affects all sectors of the dairy 
industry from the farm production of milk to international 
trade in processed dairy products. The primary goal of 
practically all dairy policy, whether under federal or 
provincial jurisdiction, has been to encourage long term 
viability in the domestic dairy industry so that domestic 
consumption of dairy products may continue to be met 
predominantly from Canadian production. The most appropriate 
avenues Of promoting this viability were perceived to be 
stable farm prices, with consequent reduction in risk, and a 
sufficient level of returns to encourage dairy enterprises 
to be efficient and adopt new cost saving technology. As D. 
Peter Stenhouse notes: 1° 

“It is doubtful whether the dairy industry would be 

such a prominent sector of the Canadian economy | 


without extensive government intervention. " 


A. The Role of the Federal Government 

At the Federal level, dairy policy translates into a 
set of regulations that control supplies and gross returns 
for industrial milk and cream used in the manufacture of 
butter, cheese and other processed dairy products. The 
policy instruments used by the federal government to achieve 


these goals include the following: 


10Stenhouse D. P. Government Policies for the Canadian Dairy 


Industry Canadian Farm Economics Vol 14 No 1-2 
February-April 1979 p1 
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1. offers to purchase processed products at minimum 
guaranteed prices so that processing plants can pass on 
some of the benefits of higher and more stable prices to 
producers, 

2. direct subsidies to producers to increase the returns 
for manufactured milk and lower the price to the 
consumer , 

3. partially subsidizing any export of surplus processed 
products such as skim milk powder, 

4. discriminatory pricing, and import tariff and non-tariff 
barriers in order to stabilize or increase prices, and 

5. quotas, penalties and holdbacks to manage the supply of 
milk coming onto the market. 

These five major areas of policy ei vement are 

administered by the Canadian Dairy Commission (CDC), an 

agency established under the Dairy Commission Act of 
1967.The Commission has implemented a complex system of 
regulations in order to achieve the first three policy 

objectives. offering to purchase dairy products at a 

stipulated price, and with manufacturers having to compete 

for milk to maintain their plant operations, these 

Quaranteed prices in effect have become minimum prices for 

all industrial milk. To arrive at the price the producer 

will receive, an assumed processors margin is subtracted 
from this guaranteed price to obtain an “estimated producers 
market return". This return coupled with a subsidy of $2.66 


per hundredweight (cwt) of milk equals the CDC target 
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support price.'' For the 1979-1980 dairy year this price 
Witt sdnouUnt LOe> tance nero ApPrh 1, -197>" this “target 
support price was indexed with the price being determined by 
a weighting of cash input prices, the consumer price index 
and a number of judgemental factors. The weightings were 
45%, 35% and 20% respectively. 13 

Policy objective 4 listed earlier is pursued by 
restricting the import of dairy products into Canada, mainly 
specialty cheese.'4 Butter is imported by the CDC to 
stabilize domestic prices which fluctuate substantially in 
response to domestic supply conditions. Since dairy products 
are on the import control list, a permit is required for any 
importation. Thus direct control of any imports is possible 
by the Federal government. Exports of dairy products are 
supported by levies assessed on the producers of industrial 
milk.'® These levies are paid into an export equalization 
Ton « 

Policy objective 5, the use of quotas, holdbacks and 
penalties to manage the supply of milk, has undergone 
substantial changes since first becoming part of the 


Comprehensive Milk Marketing Plan in i970. Under this plan, 


11 For a more extensive breakdown of payoff calculations see 
Appendix 1 

12 Canadian Dairy Commission National Dairy Program 
1979-1980 Ottawa March 1979 

13 A price increase was triggered by the index for the fall 
of 1978 but was delayed until the spring of 1979 as a result 
of the Federal government exercising its option under this 
judgemental factor. 

14 In 1978 cheese accounted for approximately 94% of dairy 
imports 

15 For the 1979-80 dairy year this will amount to $1.00 per 
cwt and 20¢ per cwt on fluid milk (class 1 milk). 
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the Canadian Milk Supply Management Committee (a committee 
composed of representatives of each of the provinces and the 
CDC) uses a program encompassing market share quota (MSQ) to 
regulate the Canadian milk supply and to provide a balance 
between domestic production and use of butterfat in 
manufactured products. These quotas are allocated to 
individual producers based on the producers production level 
during the previous year. For the coming dairy year (August 
1, 1979 to July 31, 1980) the total Canadian requirements 
for milk has been projected to be 98 million hundredweights 
of milk and with a 54 sleeve allowed to ensure this supply, 
103.7 million hundredweights of MSQ have been allotted. The 
sleeve represents production that is surplus to demand but 
needed to ensure this amount of production is forthcoming. 
To ensure against surplus supply, individual producers who 
exceed quota allottment will be subject to an over quota 
levy or penalty of $8.00. A further continguency levy of 25¢ 
per cwt will be collected on all in-quota production to 
cover any disposal costs associated with the in-sleeve or 5% 
surplus production that ultimately may prove surplus to 
Canadian requirements. It is refundable to any individual 
producer at the end of the dairy year if he has not produced 
in this sleeve or if production into the sleeve is needed to 


meet Canadian requirements. '® 


16 A more extensive review of Canadian dairy policy can be 
obtained by consulting the Stenhouse and Bray articles cited 
in the references 
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B. The Role of the Provincial Government 

Responsibility for ensuring the orderly marketing of 
milk in Alberta rests with the Alberta Dairy Control Board 
in Wetaskiwin. This Board administers the MSQ aspects of the 
Federal program, represents Alberta on the Canadian Milk 
Supply Management Committee and determines production quotas 
for the manufacture and distribution of fluid milk. Quotas 
for the fluid market are determined by the projected 
requirements of the market area. To ensure a constant 
supply, this requirement is set at approximately 115%- 
130%17 of fluid sales with daily quotas allotted to 
producers to ensure this amount of milk will be forthcoming. 
Production that is surplus to fluid requirements receives 
industrial milk payment and subsidy provided sufficient MSQ 
is held by a fluid producer to cover this surplus 
production. In recent years the amount of fluid production 
that has entered or "crossed over" eho. Un industrial milk 
market has intensified. 

In the past 4 years cross over has occurred from the 
industrial market to the fluid market as well. Industrial 
milk producers who meet certain sanitation requirements, 
daily production averages, etc.'® receive 400 pounds per day 
of fluid milk quota their first year of entry and 200 pounds 


17 Although this amount of quota is issued, payoff for any 
month is based on the level of sales, approximately 85% of 
quota allotted, with the remainder receiving Class 2 and 
Class 3 prices. 

18 For a complete list of the requirements for entry see 
Alberta Dairy Control Board Milk and Cream Quota Policies 
for the 78-79 Quota Period Wetaskiwin April 1978 p1l2-14 


: I 1 f at 7 
a - : r 
¥ = oh are 

haaowrtawe no (ersaPpyor?. ott 'S 


, i » t . 
pei 7am. vt gba) get at anit a voit tea as hs oF 
We oT 


pasod i orshod valet so tedha ef mie alge: etsentha mt AO 


| ‘ ; 02M oH) posh hatte Stead) at path Janbe 
e ‘ M Se , 7 v 
~ Lol es ie 
ait + 29S) 2Itges toe" aes TENS +i ato 
MieD. 7 
' <7 
Otc 2 ipod temapgent vi 
Hitit tm el ludtateéo brig wre ius Rayo 

ar ¢ afl 

ol 3aaTotg tet go pewirioetey 24986 rapid tent sat 
N83 ai 1s & TM < sg te TS wreds ony eine 7184 e: 

(" i 
ac7t ¥1S) oF ie ag re Set i.e eS atid tt 


a) a vf MS 


oF betlvolts esos yi tee 7ihhw asiea brine 


( 2gineme7iune> ferute’ 2? Si. rue at }ga7: net ea 


\/ it . 


i iouborg, bhult ta deems seit shay ISIE, iy i 


s'otautnd wf ont “sesve- weeeo 1g" 49 versie! pe | 
; ie 
, seh Nora ail Nii 


#1) moat? fe rtbte cat ud weed. a’ aay b 78a anit Ah 


Peru devisn bt eer #5 st Nee mn. sAd, Ke leitem LAT. i. 


s@lnewetty oe erica ra es ie bi ebm ort zracwboa 
yeh tag aibrmog O08 sien: ah ‘ote rene TONE | Oe SU ¥PF 


mbowag O98. bop. yrtia) Hy nb matt Fe Hit <a : 


\ ras 
; | . a) Rie y ee 
ne i) Bp q (beues! foe =; cad 
4 4G8 vlotant <o7age as Feiner ss tps 


sulk a eo my - Dagael ek bLikie 

ri whe 227 ¢ we) sere v tc 
‘aa 2. e — ape. ~ anna . ibdee, 
guse ho? Binexnayh pS) | 


{2 


for each of the next two years. These producers, referred to 
as graduated entrants, can now receive the higher returns of 
the fluid producers for this amount of shipment. The 
ultimate conclusion of this program appears to be one 
classification or standard for milk with milk utilization 


determining the payout to any individual dairymen. 


C. The Changes in Current Dairy Policy 

Alberta dairy policy planners have made four major 
changes to present policy to stimulate the increased 
production of milk in Alberta and there is a possible fifth 
change to come. The first major change was with respect to 
market share quota. Now 

"all existing milk and cream producers who wish to 

deliver more than their present quota can do so and 

the Board will cover such shipments with market 

share quota from the provincial reserve. This policy 

is in effect now and will be effective until March 

St d980% he 


The second major change was to the graduated entry policy. 


1. application can be made at any time (as opposed to 
applications being received prior to September 1 for 
entry January 1st of the following year), 

2. dates of entry into the program are now the ist of 


January, May and September (as opposed to only January 


198 Alberta Dairy Control Board Important Information to All 
Dairy Producers Wetaskiwin February 1979 
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3. any applicant must now have marketed for two months 
prior to the date of entry (as opposed to the previous 
12 months), a daily minimum production of 500 pounds of 
milk or milk equivalents. This amount marketed will 
determine the entrants total quota and will not exceed 
800 ibs, 

4. all applicants must be able to show that they have not 
disposed of either fluid quota or quota previously 
allotted under the graduated entry program during the 24 
month period prior to entry, and 

5. all applicants now are required to ship to their 
graduated entry quota for a full two years after final 
allottment before such quota may be transferred other 
than as set out in other regulations. 2° 

A third change and one yet to take effect is the advising of 

fluid quota holders that commencing in 1980, September is to 

be considered as a permanent low month. Previously, the 
three lowest months of the year when the ratio of total 
provincial production to fluid sales was lowest, was used to 
adjust fluid milk quotas. 

The fourth major change to policy has taken place in the 

area of the Alberta Dairy Development Program.?' A fifth 

major stimulus to the industry could be the reintroduction 
of the Dairy Heifer Import program, a program that was in 


place for 2 years in the mid 1970s. 


OVO Diao . 
21 See Appendix 2 for complete details of this program. 
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D. The Economic Theory of Milk Supply 

The rationale behind this study and the assumptions 
implied in making the dairy policy changes listed are all 
based on underlying economic concepts. In a mixed economy 
such as Canada’s, the prices paid or received for goods and 
services are set by the forces of supply and demand and 
therefore subject to the choices and desires of producers 
and consumers. Within the dairy industry, a free market 
economy does not totally exist as prices that producers 
receive for their milk and the prices consumers pay for the 
product are controlled by Public Utility Boards. In 
addition, milK quotas act as as barrier to entry into the 
industry and restrict the supply of milk coming onto the 
market. 

In Alberta these quotas form only one causal variable 
in determining the aggregate supply of milk. This aggregate 
supply function can be defined as the relationship between 
the quantity of milk supplied and those variables that 
influence that supply. Symbolically such a relationship can 


be viewed as follows: 


where: 
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Q = the quantity of milk supplied to the market 
per time period 
C = the number of cows-in-milk per time period 


T = the technology employed in producing milk 


P = the current price for milk 

S = the prices of other products 

I = the prices of imports 

O = the objectives of dairymen 

N = the number of dairy farms in Alberta 

R = the size distribution of dairy farmers supplying 


the market place 
G = policy variables (support programs, deficiency 


payments, quotas, imports, etc.) 


In this study the supply function assumes the latter 8 
variables to be independent of quantity and any movement of 
these variables to be stochastic or random in nature. Milk 
production is therefore theorized to be a function of the 
number of cows-in-milk in the Province. This relationship is 
assumed to hold only in the short run since a long run 
analysis would require the inclusion of the other variables 
mentioned above. Since no accurate data exists on the total 
number of cows being milked, a sampling of herds is used to 
estimate the average production for a herd and then this 
estimated production is multiplied by the number of herds in 
the Province. 


At the individual farm level this average production is 
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obtained by means of a production function with the number 
of cows-in-milk as the major input variable. In the ideal 
economic setting this can be explained as follows. Dairymen 
combine their major input factors (land, labour, capital and 
management) to produce milk. Dairy cows as a factor of 
production form a part of the much larger capital input. If 
it is assumed that all inputs are fixed with only the number 
of cows free to vary for a given dairy farm, a relationship 
exists where milk yield increases as more cows are freshened 
or added to the milking base. As more and more cows are 
added, the law of diminishing returns eventually becomes 
prominent. This law states that as “increasing amounts of 
one input are added to a production process while all other 
inputs are held constant, the amount of output added per 
unit of variable input will eventually decrease." As more 
cows are added to the herd, adverse effects are felt from 
other constraints such as management skills, technology, 
etc. At such a time, the marginal product of cows or the 
increase in total milk yield with an addition of another cow 
will be zero. Correspondingly the value of its marginal 
product will be zero. 

This point is the upper boundary at which cows are 
added to the herd. However before this point is reached, the 
dairymen would have reacted by altering al! or some of the 
input constraints that were previously assumed constant. The 
least number of cows which ought to be employed in the herd 


is that number which results in the highest average product 
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produced per cow. Thus this same dairyman would keep adding 
cows to his milking base if the additional product that each 
additional cow contributed was increasing. When the lower 
boundary point is reached, the additional product each 
additional cow contributed equals the average product of all 
the cows. The area bounded by this upper and lower boundary 
points is referred to as the Stage Two area of the 
production function and is the only rational area of 
production. 

At the individual dairy cow level, there is an 
under lying production function unique to each cow. This 
uniqueness is due to the inherent differences that exist in 
the productive ability of cows, nutrient requirements, 
stomach capacities, etc. If two inputs such as grain and hay 
are viewed, and all other inputs are held constant, a locus 
of all the combinations of grain and hay could be shown by 
means of a product curve or isoquant. The tradeoff or rate 
of technical substitution at which grain is substituted for 
hay to maintain a given output level is determined by the 
price ratio between these two inputs. More grain would be 
substituted for hay if the price of hay was to rise. 

There is also a tradeoff between cows-in-milk and 
cows-dry at the farm level and this tradeoff coupled with 
the tradeoff between production inputs for each cow 
describes the fundamental decisions that face an individual 
dairyman. The dairymen’s optimizing behaviour is to vary the 


levels of both cost and output so that profit is maximized. 


aitiohs gear btudw migihyanek: dina an bala. is “won 
jauisorg bene? eGbe set. 7% eR ene : int cat 
Latent oe panty gn hesetae Qt ey aye ee: t a 
toubor Tei 2roee ers beorlngied: at. twee 
ibere), aps 7sys et 2heupe betucat tray WOO, 301 
oun! Ue TeGd ani: wa westorniad sorts weit on 


1s OWT. opet2 ort ob a) Speiste af, 


<& 


s2>o6 fartoried vital. ef one not fomu? net: 

ma @7 ateas .pensl wes ery nies of Pepe art Aan 
sit? .who eee ot sus ttul nat oat aoe ort ed: , 

nt 2@9xe teas eaometaly'b (nbseanh ere oh, aut et | . 
| . 
pe Ionstaiun awe to nee ieee uaa . 
pata 
“or OAs nresD, a6 fous: sucNtt. mie oe), eee (260s LoBaD 
aiisal 4 insletiog Gian. e976" e) Jah dria ite Orie. oad 
d pwieta el blues yet bbs nl Bag tis Prot T ph cele ort ih 
1 RESO) aK T og Sper E oe san Powho wl a — 
feds si idherp tape fe wiot Aus, edu Vea 
ent ye inteigb, of Deved -tigeige nr yg § biserneculs: 
20 Divow 7 stgHgmOl 8 huseint RY aarp mew ted viene 
eae tot eae «eA to5aae He *) Yee: sat beaté 


so ~ << © 
O89 4 ur 


a 


beg 4) tee ~embo regent hoes b! 901s eh reat 

Hbtw abiquas ++dabree?! ius hey pddet tess? anit’ deg 
van fisge to? etuanh sobisubanay notheted oak 
heonivihe! me epee heed atwisiony eal etd. ig 
add yay of 2t. Quntyeced .onisl 


{boride Aen ef .t2team tedd ga: 


18 


Total revenue is determined by multiplying the market price 
of milk by the number of hundredweights of milk sold. Costs 
are computed on the basis of equating the per unit cost of 
each input in the herd with the value that each input adds 
to marginal product. 

The approach that Alberta dairy policy planners have 
taken to stimulate the production of milk draws heavily on 
this optimizing behaviour. The change in market share quota 
whereby any overproduction will be covered by quota issued 
by the Alberta Dairy Control Board is an attempt to remove a 
supply constraint that may have forced some dairy producers 
to produce at a point where the marginal factor cost was not 
equal to the value of marginal product. 

In an output sense this constraint can be shown as a 
point where marginal revenue (MR) exceeds marginal cost 
(MC). Diagrammatically this relationship is shown in Figure 
2.1. At price (P) the dairy producer would like to produce 
Q*, the point where MR equals MC. However quota restrictions 
allows only the quantity Q1 to be delivered. The area shown 
as ABCD (the area bounded between the average cost of 
production and the price of output for quantity Q1) 
therefore represents the excess profits that accrue to the 
producer as a result of the enforced underproduction. A 
result of this action is the capitalization of these future 
anticipated profits into the value of quota. Removal of the 
quota restraint should in the long term allow Q* to be 


attained with its attendant increased milk supply for 
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Figure 2.1 Milk Production under a Quota Constraint 


S/Tb.) 
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Quantity of Milk (1b.) 
Albertans. 

The easing of requirements for graduated entrant 
admission allows for increased revenues currently only 
enjoyed by the fluid milk producer. These higher returns may 
also induce more producers to enter the industry or convert 
from cream shippers to milk shippers. The changes made to 
the Dairy Development Program is also an attempt to lower 
some of the high costs associated with the industry. 

This brief overview of the relevant economic theory not 
only forms the basis for this study and the recent policy 
changes but also can be discerned when previous work done on 
predicting milk is reviewed. This previous work will now be 


discussed. 
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Peer reRATURE REVIEW 


A. Previous Studies Related to Predicting Milk Production 

A number of studies have been undertaken in the area of 
predicting milk production. They can broadly be classified 
into two categories; those predicting milk production for an 
individual cow for use as a farm management tool and those 
attemping to construct larger more complex models to 
forecast aggregate milk production (supply response 
analysis). Each of these study areas have used different 


approaches to predicting milk production. 


B. Predicting Milk Production for Individual Cows 

Work involved in predicting milk production for an 
individual cow and, when aggregated, for a herd has revolved 
around demonstrating how dairymen could make better 
management decisions by using partial days in milk records 
that were extended or extrapolated to 305 days. The most 
important aspect of predicting at this level has been in 
estimating the shape of the lactation curve for a particular 
cow in milk. Extensive work in this area has occurred since 
the early work of Sanders in 1930. After analyzing nearly 
4000 lactations for cows of various breeds and ages he 
concluded that the curve could be regarded as a combination 
of four elements: total yield, maximum yield, persistency 
(or production falloff through time) and spring hump 


seasonality. Much of the work done to date has substantiated 


20 


OSV! GvVva" wari io) 5 Is? waht pring nae ors: WD rounii 


/ ‘ Pr ive 7* ‘ 7 
‘sid asgn Ons abae rd ote esos mye ’ 
_ Seb tent dma & zerbebaspen - ote ev 1D eit Yar iz 


Oy pews ah it ohahy naa a ibt 


WITVah SROTAMa te. PP 


| ae Mm ti —~ 
' ve . . " i AY A 7. 

nt dowhond EN ont IRN oi bots hai an haya? auol he . 
. A. 

i a? (46) °SRAL) (esi avid eertite oO ‘vee min A ‘¢ 

is ~ i) satsah ae eras 

59 Vert! . cee > Sow ai hate ani toi 


: + seortl j2e% nce tas out 
'> ~t : a phic. 
es SBS Ge (Ss? € 26 @eh i wd — bh 


J 


Ye! MOD SA hate | yk 3 avai oy orien nied 
7 ¢ 
i sf , ‘ 


lia) mot yom a. a) Toes isla ims | jenoen I 
Ae TeTtib besu eved seach ybia eeenl 6. fond ted 
Te) ToJaote A an Ot afte ef arse 


1 oh Lg ‘f it : 


suite? ic Was i aan Ee tart igibent iy 
GF Pouhowws Ml tie Lf ae3 3 ym ne bavi ownt ariel = =e 


PRG a 


‘Gen ste 


vetted. sManm bhups cyanea? eo wel Set ett ericomota 
torXacey A hat At eyed Part ney pil ats e ‘tei a toe: an re - 
igom ant  -exah 6 ge ashing ” hoon tase 


| F } ha . | 
ny need 22 ewer oi Hy Diy * #0 te) q - we bins ais 
| oe ree 
Si Yort160a-5 ef-s av abe rt rat; 


. to ae ee i 
sandra baa Tud9e6 alah asort anes at he ovr enate? a 


i see haleh Ras. 1610) ORGT see. ej 


t on ited ; ru a 


oi. anek 
. 


» 


| - pth ae Byte ve my 


rat 


21 


these findings with major work done in showing the shape of 
the curve and the accuracy and scope with which it can be 
used. Initial attempts to use the characteristic shape of 
the lactation curve for predictive purposes involved 
determining ratio factors to extrapolate 305 day production 
from part lactations. Lamb and McGilliard (1960). developed 
such factors to investigate the relative effects of breed, 
herd, age, parity?2 season of freshening on the relationship 
between partial and complete lactations. Their findings 
suggested a difference existed between heifers and older 
cows. Although this finding was substantiated by later 
research it was not in complete agreement with the 
literature of the time. A year later Van Vleck and Henderson 
used regression analysis to extend partial peeteah. 13 305 
day completion. Their main interest was to show that the 
accuracy of predictive factors could be improved if 
intra-herd effects were used in their development although 
they concluded the extra computational difficulties 
encountered precluded their use. Other work by Appleman, 
Musgrave and Morrison (1967) supported these findings. The 
use of days open and days dry as important influences on the. 
milk production of individual Holsteins was done by 
Schaeffer and Henderson (1972). They found the relationships 
accounting for highly variable predictive results difficult 
to assess due to problems in determining if cows produced 


poorly due to short open periods or if low producing cows 


22 Parity represents the lactation number for a particular 
Cow 
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received short open periods because of their low milk yield. 
They also concluded that the effect of days dry and days 
open on milk production was largely environmental. As a 
result, adjustments of milk records for days dry was not 
recommended while adjustment of 305 day records for days 
open appeared necessary. Applemen et al showed that factors 
developed by regression analysis were more precise and 
methods were now being sought to develop more accurate 
predictive qualities by using the last available test period 
to project the remaining portion of a cow's lactation. This 
work was done by Miller in 1973. By 1973 extension factors 
had progressed to the point where they differed from those 
currently in use by the United States Department of 
Agriculture and developed by McDaniel et al (1965). These 
latter factors had not taken into account season of 
freshening or age within lactation and had also been 
developed solely from regional data. Keown and Van Vleck 
(1973) developed the most involved piece of research work 
into predicing milk production for an individual cow. They 
initally took a lactation curve and divided it into 19 
segments. Superimposing three season groupings and two ages 
for each of the first three lactions and a single age for 
all lactations of four years and over, they were able to 
calculate average daily milk yield for 114 season-age-stages 
using multiple regression. Using this procedure they 
depicted the lactation curve as having three distinct parts. 


Part 1 showed a production pattern increasing as a 


aengte-spe- ee Aft swt ftety Mitel int, eine 


re A ‘Si 
at re ae > 


iiety Alim, wat sien to soured ghar ~~ 
bas yoo eveh. ta Pomtin mele badd 
wa . cnn ve Yhap tel) enw voto a 
Ton abw Yih 246 «oF 2Zanane4 i ia 1 oF | 
cub wy eosqees wee Goce 2G d evcivid cabin oF ht 
teri? bdwede) fs Je a N dncoeerier 

s #apeaig son ey, efay leas noteseeen war 

sip iwols 910m Gohevan et, Fapwos a Sr WOR - o1gH) e 


sen a 


nae ” we 


feet sidelTavs tagt grt naahe a aoe i Lan eves 

4 Nottelos! @ ‘woo Bey Pecan Dev Ort genes om} 3 “so 
’ 

“@jee £78 yb) MTRT ee ee aa ah Ae 
no S » ved? ocechw Taha ae et Seas | 


fo Tem) dened esiel2 agHtnl ii ome ak she 
seen (888) fs To Jabadtg’ 2a Beaaiieyeh bem erul fe 
| 0s walod veget, fom-bart ecotast 8 

naa dats Hed has. abttigdseh: hit iw wpe ae ; ; 
Hoot oie 16 woe sek. Tene ppe HOP. gow | me | 
iow Notesasn th, toni Weylaver P dicge ay basanuni te 

Yad? .woo lawbFytent ae 4é¥eh@Ptagltien at Fav: gente nit 

Sit St Babi eae aye) code lant a ‘iia cia 
Zaps -owlt bane eq" mdr} hide obey naa blaga 
10% ahs, shgte othe erin! age sends jerk? ei — 
Of 34 OS S16 eB”, tee bts emg que Foam 


wena leitpisaiie at At gonial Wokasengen © 
2h 9c, forth tet peat! onived eagle 1 net? stom 
s as & gilasorany ~~? bd eatin 8 


~ : — 


as 


decreasing rate up to peak production which occurred between 
50 - 60 days in milk. Part 2 showed production decreasing at 
a nearly constant rate once the peak had been passed. Part 3 
displayed a shape to the curve that decreased at an 
increasing rate. Figure 3.1 reflects these findings. They 
also found that cows that calved earlier did not peak as 
high as older cows calving for the same season. The season 
factor was also found to be important as cows that calved in 
the first four months of the year reached a higher peak in 
production. The last four months of the year were found to 
result in the lowest peak production. These three distinct 
areas of the lactation curve substantiated the work of Wood 
(1969) and the work he had conducted as a member of the 
British Milk Marketing Board in Surrey, England. His earlier 
work had included the fitting of a logarithmic mathematical 
function to the shape of the lactation curve. 

Perhaps the most questionable work done in estimating 
the shape of lactation curves was that done by Kellogg, 
Urquhart, and Ortega (1977). Using Wood’s mathematical model 
for expressing the lactation curve they used regression 
analysis to test the goodness of fit with data obtained from 
73 Holsteins in each of the first four lactations. However 
they tested the goodness of fit after they had averaged the 
daily milk yield for each lactation. Not surprisingly they 
found the fit to be extremely good. Other work in 1977 by 
Schaeffer, Minder, McMillan and Burnside, and by Wiggans and 


Van Vleck (1979) were aimed at improving milk predicting 
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Figure 3.1 The Shape of the Lactation Curve 


19 
{ 
a | 
a7 
~ 
jer | 
Oe ; 
© 
a ’ 
13 
| 
| | 
11 
} | 
| 
| | 
| | 
PART1 | PART 2 IPART 3 
| 
7 | 


Baie She L Bee ee 


elie 2 Be a Oe ol Cee 
20 40 60 80 100 120 140 160 180 200 220 240 260 280 300315 
DAYS IN MILK 


qualities by using nonlinear models. In the latter case, 
this work involved using the last test day sample and using 
separate mathematical relationships to describe the three 
major parts of the lactation curve. 

The discussion to this point has centered around 
constructing more accurate milk predicting models in the 
face of a highly variable lactation curve. However 
extensions of a model’s ability to predict milk production 


and serve as an aid in making farm management decisions has 
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also been done. Selection, culling and mating decisions have 
been based on cows part records being extrapolated to full 
lactation. Oltenacu and Ainslie (1976), using the milk 
predicting model of Keown and Van Vieck (1973) and 
forecasted milk prices, were able to demonstrate that 
projecting the income from milk sales allowed dairymen to 
better plan their cash flow. Wiggans and Van Vleck (1978) 
demonstrated that incorporating the proportion of 
concentrates fed to dairy cattle into a milk predicting 
model did not alter the shape of the lactation curve and 
hence its use in predicting milk for farm management 
purposes. 

The most extensive use of the shape of the lactation 
curve for use as a farm aid is the Herd Management Control 
program administered by the British Milk Marketing Board. 
This program makes use of the vast amount of information 
held in the computer based National Milk Records and is used 
as a management aid for consulting officers and dairy 
farmers. As R. J. Clothier’s summary of the program to the 
Edinburgh Dairy Conference in 1977 stated: 

"By comparing each cow’s predicted yield with her 
recorded yield at each monthly recording it is 
possible to identify the problem cows, isolate the 
causes of poor production and plan corrective 
action. The sum of individual predictions provides a 
reliable prediction of herd yield. 


This herd yield is plotted on a graph. The 
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graph serves three main functions. Firstly, it acts 
as a check on herd production from week to week 
between milk recordings as the farmer is expected to 
plot weekly yield against predicted yield and to 
note the relative slopes of the two lines. Secondly, 
it provides a picture of the years production 
against predicted production which aids in 
identifying periods when management could be 
improved. Thirdly, the graph may serve to monitor 
changes in management. " 

The basis for this service has been developed from the 
work of Wood (1969) and his analysis of the lactation curve. 
In this program the management service serves not so much as 
an accurate predictor of milk production but as a management 
aid in alerting the dairymen of the need to reassess his 
management practises, feeding schedules and quality, etc. if 
the graph slopes between actual and predicted milk yield are 


not the same. 


C. Predicting Milk on an (eerecne Basis 

Studies in this area have aimed at evaluating industry 
supply in response to an increase in milk prices since 
dairying is a highly regulated industry. The hypothesis used 
was that estimating this response accurately would provide 
dairy policy makers with the ability to provide for an 
adequate supply of milk and help maintain a "reasonable" 


balance between production and consumption. Brandow (1951) 
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used simple single equation regression analysis to estimate 
milk supply for the United States. Halverson (1955) using 
the same statistical procedure expanded this work further by 
estimating milk production per cow as a function of the 
milk-feed price ratio, hay production, and cow numbers. This 
work was conducted on each of six regions and in all cases 
production was found to be highly inelastic. Halvorson 
(1958) used a Nerlovian distributed lag model to determine 
the short and long run price elasticities for milk in the 
United States. He found these to be in the ranges of 0.15 to 
0.30 and 0.35 to 0.50 respectively. Ladd and Winter (1961), 
using single equation regression analysis on time series 
data, estimated milk supply elasticities as a function of 
average output per cow and the number of cows in the United 
States on January 1 of the year of prediction. Other 
variables were encompassed into the model to properly expand 
upon the estimation of these two major variables. 

Wilson and Thompson (1967) used simultaneous equation 
models to determine equilibrium price when both supply and 
demand relationships were considered. Zepp and McAlexander 
(1969) reported that regression analysis was more 
appropriate to use than non-flexible price lag structures. 
Chen, Courtenay and Schmitz (1972) using a polynomial lag 
formulation and the premise that the quantity of milk 
produced in a given time period is a function of the price 
of milk, prices of inputs used in milk production, returns 


obtainable from competing commodities and the existing 
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technology, measured the lagged output response resulting 
from a changed product price. They concluded that output 
response to a given price change increased through time and 
then declined; that a polynomial lag formulation offered 
added advantages over that of geometrically declining 
distributed lag price structures; and that in general the 
whole area of supply response was fraught with certain 
estimating and theoretical problems. 

To date the use of distributed lag models to predict 
milk production appears to be a research area open to 
different interpretations. Firstly, Griliches has stated 
that most distributed lag models have almost no or only a 
very weak theoretical underpinning. Thus it is hard to 
determing which lag models to accept or reject. Secondly, 
the use of time series data in which more than one price 
change may have occurred results in measuring the 
responsiveness of production that is not predicated on a 
once-and-for all price change. In view of these limitations, 
measuring milk production response using distributed lag 
models has generally been viewed as less than a perfect 


science. 


D. Relevance of the Literature to this Study 

In attempting to predict milk production in Alberta, 
both the work done on the shape and characteristics of the 
lactation curve and the concepts underlying supply response 


were drawn upon. The development of the four models for this 
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study had to meet two criteria; 1) ease of use for 
predicting supply and 2) use the Dairy Herd Improvement 
program as the basis of study. Alberta milk supply was 
theorized to be a function of the number of cows-in-milk for 
a given month. This simple functional relationship was 
expanded upon in the Complex models by segregating dairy 
herds into common producer categories. The work of the 
lactation curve was used in the development of the 
Stage-of-Lactation model to determine if the number of 
cows-in-milk at a common part of the lactation curve for a 
particular producer category would improve on the 
predictability involved as well as serve as a ready 
"barometer" of cow movements. To the extent that aggregate 
milk production is only a function of the number of 
cows-in-milk and that other variables that might suitably be 
appropriate for measuring supply response are viewed as 
being stochastic, the overall accuracy of the prediction 
might be hurt. However the absence of these variables allows 
for ease of use and may prove to be quite acceptable if 


accuracy of use is not affected substantially. 
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IV. METHODS AND PROCEDURE 

The main objective of this study was to develop an 
econometric model suitable for predicting milk production in 
the Province and to compare them for accuracy and ease of 
use for the eight month period November, 1978 to June, 1979. 
Four models were chosen for development to reflect two 
sample bases and two sampling procedures that were assumed 
to the study. For ease of use these four models were 
identified as Complex and Simplified Stage-ofLactation and 
Complex and Simplified Cows-in-Milk models. The only 
difference between the complex and simplified designations 
was with respect to the sample size and the designation of 
individual producers by type of shipment. The simplified 
models dispensed with this latter differentiation and were 
developed from a smaller sample base. 

The essential element to predicting milk production for 
a particular month was to predict average monthly production 
for a herd on DHI and then extrapolate this average 
production over all DHI herds that were of the same producer 
category. To encompass those herds that were in the same 
category but not a participant of the DHI program, 
deliveries between DHI and non DHI herds per month would be 
compared over the 12 month period of July 1978 to June 1979 
to obtain a monthly production coefficient. Total production 
for a particular category would be the sum of the predicted 
production of DHI herds plus predicted production of non DHI 


herds. This aggregate production would then be multiplied by 
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a monthly delivery coefficient which represented the 
proportion of milk produced in a given category that was 
actually delivered. Summing over all five categories 
determined total provincial milk delivered for a given 
month. 

To select the production categories seasonal production 
patterns for different types of producers were observed over 
the past three years. Although milk producers could be 
categorized into simply fluid or industrial producers, the 
cross over between the two markets in the last few years 
suggested more well defined categories might be advantageous 
in order to achieve increased accuracy. It was theorized 
that fluid production, because of its daily quotas, should 
have an even seasonal pattern of production while industrial 
producers were theorized to have a curvilinear pattern of 
production with production peaking in the May-July period 
and Bete lowest in the November-December period. The two 
classes of graduated entrant were theorized to behave 
somewhere between these two extremes. The Graduated Entrant 
1975-77 approximates a fluid producer due to the 800 pounds 
of daily milk quota earned as a result of being in the plan. 
In addition, additional fluid quota can be earned if 
production takes place above quota and more fluid quota is 
distributed as a result of increased provincial consumption. 
The Graduated Entrant 1977-79 is locked-in with respect to 
fluid milk quota and any overproduction enters the 


industrial milk market. Hence this type of producer was 
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theorized to behave as an industrial producer. Therefore 
each milk producer could be classified into one of five 
categories. These were designated as Traditional Fluid, 
Graduated Entry 1975-1977, Graduated Entry 1978-1979, small 
Industrial milk (less than 10200 Ibs MSQ) and large 
Industrial milk (greater than 10500 Ibs MSQ). To properly 
identify each producer’s category, a computer printout of 
all dairy producers in Alberta was obtained for June 1979. 
Each producer was designated into one of eight different 
categories as well as any participation in the Alberta DHI 
program. Because not all producers on the DHI program could 
be properly identified from this information, a cross check 
of all producers in the DHI Herd Master file was made with 
the Alberta Dairy Control Board’s Market Share Quota List. 
Cream producers were disregarded for purposes of this study 
and the two small industrial producer categories were 
amalgamated into a single category. A sampling of DHI and 
non DHI herds from each category was done for development of 
the respective models and for comparison of milk deliveries 
within category. The number of producers in the Province by 
category, the number of DHI herds sampled for the study and 
the number of non DHI herds sampled for calculation of the 
monthly production coefficients is shown in Table 4.1. 

The sample size was determined by calculating the 
standard deviation of production for each category chosen in 
the study and selecting a standard error that was acceptable 


for the study at a .05 significance level. Standard 
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TABLE 4.1 


JUNE 1979°DAIRY PRODUCERS BY CATEGORY 


Category Total Total Samp led Samp led 
in Prov. on DHI (DHI ) (non DHI) 

Fluid 636 283 62 30 

GEI1975-77 157 (eS) 25 30 

GE 1978-79 302 101 26 30 


Industrial 


less 6500 210 le 7 1S) 
6500-10500 215 Ag: 8 15 
more 10500 188 47 19 30 
Cream 

less 500 1099 2 2 - 
500-6500 2874 19 a = 
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deviations for each category of producers was calculated 
from an Alberta Agriculture computer printout for DHI herds 
in each category. The printout contained data for the period 
November 1977 to October 1978. Although this printout was 
incomplete, it did estimate the mean and standard deviation 
for each category. The mean and standard deviation for each 
category for the month of October is shown in Table 4.2. The 
high standard deviations for each category relative to the 
means suggested a large sample size for each category should 
be used if the .05 significance level was desired. Hence 
assuming a normal population, a total allowable standard 
error of 4000 kgs., a significance level of .05, and the 
October herd printout of standard deviations the following 
formula was used to calculate the sample size for each 


category. 


AQ001/*t296eane <ScDem/ ssquaré mootsof en 


4000 = the allowable standard error for each category 
1.96 = the confidence interval of 2 standard deviates z 
S.D. = the sample standard deviation of each category 


N = the number of herds to be sampled 


This large sample size was possible for only the fluid milk 
category since the other four categories had too few herds 
enrolled on the DHI program to meet the determined sample 


size. The number of herds calculated for the sample was 
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TABLE 4.2 


DHI OCTOBER PRODUCTION BY CATEGORY 


se S- — —_- 


Category N Mean Standard Deviation 
Fluid VAT Si3/7 17602 
Gees SHOR ss 76 eov2en 8888 
GE 1978 40 20289 weLg 
Industrial 107 ARS 9030 
Cream 29 7491 6449 


subject to further change as a result of the inability to 
cross index DHI herd numbers with the corresponding numbers 
of the Dairy Control Board printout. 

To check the production assumptions regarding the 
producer classifications chosen as well as assist in the 
selection of an appropriate model specification to use, 
average monthly production for the sampled herds in each 
category was compared. Figure 4.1 shows the results of these 
comparisons since December 1975. 

Since the DHI program is a computer based program, 
tapes of monthly production data for sone sees TES on the 
program for the time period November 1978 to June 1979 were 
obtained. For each month this consisted of herd statistics 


(a Herd Master file) and individual cow data (a Cow Master 
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To develop models that considered different aspects of 
cows-in-milk, both the complex and simplified models were 
developed from each data file. These were respectively 
called the Cows-in-Milk and Stage-of-Lactation models. It 
was theorized that the Stage-of-Lactation model with its 
more extensive breakdown of cows-in-milk into six different 
lactation stages should be of greater benefit to policy 
planners. The simplified models were developed from a 
sampling modification of that used for the previous models. 
They were developed to see if the increased ease of use was 


more than offset by the expected reduced accuracy. 


A. Testing For An Appropriate Model 

It was theorized that a linear model specification 
should be the most appropriate model to use since the 
relationship between a herd’s production and the aggregate 
number of cows in the herd should be linear in nature. To | 
test this hypothesis the Stage-of-Lactation model was used 
with the Traditional Fluid and small Industrial-milk 
categories (these two categories represented the two 
production extremes shown in Figure 4.1). Linear, double 
log, and two semi-log models were tested with November’ s 
stages of lactation predicting milk production for each 
month up to June. Selection of the best model to use was 
predicated on the coefficient of determination (R2), and the 


standard error of the estimate 


+> wWleeae: Instr he welanae vert efe : 
ae An a ; 

aa tf oor ror ihe. el wrod, arts woh iw ns 

7 7 ag fais a 

iavitseqees svaw sige 6th 21s) cone ay Bee 

- piebon Norieiasi-Feragere vee. ay iret “ee. 
itiiw been mot yey 26+ to> eget ans “tai ete a 
wan fa aE ATT eM 7 nwooNnes td sie 
Texse to ed bl uote, sepela nate 
wor et | avah a sw af ebbcen net? ragits a ie" * on 

: mm CoO TvS ty «Waa fat) saris ba) norinstitteneant f n 
a8w dati to gase heggeniods el 4) see bt neqe! svab -a-yaMe \ at 
| . a4 
(oswol8 benuhe 4. baioeaae O99 yo Joetto nent 


ad ay 


i 


Od Ve or ei per) 


} 


fsbo% ani honergya vw 704 wien 
noi leat ?igege lebam nseenft i« Yate ektsoer! esw HL: 
wm? 4 we S64) oO} Tebum idl i jaon at ad ‘itis 


- ,: 
MW 


stenscnge acld bee Act toubore ae ee 8, Hemited phe ot 
of, .s1gree ny teenit nis uf Liga Gabe! On. ve eos Kind 
bgeu asw ($nom! Hol siagd + a anh * eatt ubverttedat ary? . 
Ait *tauhri i tie ive Whur? iusto a 
- oe ont ad Tee Fe Faihtgeres: tw? qapnt) its | 160, 
Sfdueb. .aeanh) Whe nde id aac enwentae hh b "g 4 

& “Weamavan sy tv bo tad hie ai stig pot Tene aan of 
Hoaes not moifaubae oH tm onhtér eer noi satan We a a 
aim on Gt Ighom 1ée¢ att) Yo not last ee ae af guns 

ont ons .(*h) notseninmeteh Fo inetatttecs self 


Sattss 94? to 40 


38 


B. The Complex Stage-of-Lactation Model 
Data for each cow in a particular herd is stored in the 
DHI Cow Master file. This data includes such information as 
her age, number of days in milk, days dry, days in calf, 
etc.*3This model was developed using the number of days in 
milk as the important determinant of production for each 
category. It was theorized that using different stages of 
lactation should be more descriptive of herd changes and 
supply the user with more information when actual and 
predicted yields differed significantly. Drawing upon the 
literature, the lactation period of 305 days was broken down 
into 6 different stages (dried off was included as a stage). 
These were chosen for the following reasons: 
0-40 days - this stage would reflect additional 
-calvings immediately 
41-80 days - this stage would reflect the peak 
production of cows 
81-150 days- this stage would reflect any falloff 
after peak production was reached and 
still at a high level 
151-305 days- this stage would reflect on cows 
declining in production and should be 
linear in nature 
306+ days- this stage will reflect cows 
falling outside the normal lactation 


period and having minimal production 


23 See Appendix 3 for a printout of the complete information 
available on each cow. 
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Cows Dry -_ this stage would be viewed as a "slack" 
variable and should reflect possible milk 
production response three or more months 
into the future 

After sampling from each category in accordance with 
Table 4.1, data from the cow master file for all the sampled 
herds was separated into each of the six stages listed 
above. Ordinary least squares (multiple regression) analysis 
was then used to derive prediction equations for the same 
month and from one to as many seven months into the future. 
In its linear form the model for an individual herd in a 


specified category would appear as follows: 


YOjim) = "BO feb 1 Cm) AN ee eb 2 tm X27) +7 BS tmyeRs (a) + 
B4(m)*X4(j) + B5(m)*X5(j) + B6(m)*X6(j) 


where: 

Y = monthly production in the jtmth month where 
j refers to the data or base month and m refers 
to the number of months Reine predicted into the 
future 

BO = a constant 

B 1 = a coefficient associated with the number of cows 
in milk for the 0-40 day period 

X1 = the number of cows-in-milk for the 1 - 40 day 


period in the jth month 
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B2 = a coefficient associated with the number of cows 
in milk for the 41-80 day period 

X2 = the number of cows-in-milk for the 41 - 80 day 
period in the jth month 

B3 = a coefficient associated with the number of cows 
in milk for the 81-150 day period 

X3 = the number of cows-in-milk for the 81 - 150 day 
period in the jth month 

B4 = a coefficient associated with the number of cows 
in milk for the 151-305 day period 

X4 = the number of cows-in-milk for the 151 - 305 day 
period in the jth month 

B5 = a coefficient associated with the number of cows 
in milk for greater than 305 days 

X5 = the number of cows-in-milk for the 306+ 
day period in the jth month 

B6 = a coefficient associated with the number of dry 
cows 

X6 = the number of dry cows in the jth month 


This process was repeated for each of the 5 categories 
until five equations were obtained for each month being 
predicted. Coupled with 8 equations for each category for 
November, 7 equations for December,etc.,........ and 4. Tor 
June, a total of 180 equations were obtained for this model 


over the time period under study. 
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The equations predicted average production for a DHI 
herd for a given month. A monthly table was constructed from 
a comparison of milk deliveries to include those producers 
outside the DHI program but within the same category. This 
table was produced by randomly sampling herds outside the 
DHI program and comparing milk deliveries with those DHI 
herds used in the model development. This was done by 
determining the dairy plant to which each sampled producer 
shipped. Comparing deliveries over a 12 month period and 
calculating the average monthly shipment for each category 
allowed determination of a monthly production coefficient. 

Therefore the total milk produced for a specific 
category in the Complex Stage-of-Lactation model could be 
expressed as foie 


Y(jtm) = N*P(j,K) + Gxen*P( j,k) 


where: 

Y = the total milk production for the kth category in 
the j+mth month where j = the base month, m = 0,1, 
Pe Tari g ol er tad ee ere 5 

P = the average predicted production of herds on the 
DHI program for the kth category and jth month 

G = a production coefficient for dairy herds outside 
the DHI program in the mth month 

N = the number of herds in the kth category and jth 
month enrolled on the DHI program 

n = the number of herds in the kth category and jth 


month not enrolled on the DHI program 
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The amount of this milk that was actually delivered to 
plants could be calculated by multiplying total predicted 
production for each category by an appropriate monthly 
delivery coefficent. Total predicted deliveries for a future 
month could be determined by summing this predicted milk 


delivered over all the categories for a given month. 


C. The Complex Cows-in-Milk Model 

The only difference between the Stage-of-Lactation 
model and the Cows-in-Milk model was in the degree of 
complexity with which the number of cows-in-milk was 
expressed. The data in the herd master file consisted of 
each herd’s production, number of cows in the herd, and 
percentage of cows-in-milk by month. Appendix 4 shows the 
typical information in this file. For this model, monthly 
production was theorized to be a function of two variables: 
the number of dry and milking cows in the herd at any given 
time. In its simple linear form this relationship for a 


particular herd in a specified category would appear as 


follows: 
Y(j+m) = BO + B1(m)*X1(j) + B2(m)*xX2(j) 
where: 
‘ = monthly production in the j+mth month where 


j refers to the data or base month and m refers 
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to the number of months being predicted into the 
future 
a constant 
a coefficient associated with the number of cows 
in milk 
the number of cows-in-milk in the jth month and 
Kth category 
a coefficient associated with the number of cows 
in milk in the jth month and kth category 
the number of cows in the ith herd in the jth 


month 


The overall Herd Master Model for predicting milk production 


in the Province could therefore be expressed as follows: 


where: 


Y(j+m) = N*P(j,K) + Gxn*P( j,k) 


the total milk production for the kth category in 
the j+mth month where j = the base month, m = 0,1, 
2, POM Verand ak 5 Sipir ot: 5 

the average predicted production of herds on the 
DHI program for the kth category and jth month 

a production coefficient for dairy herds outside 
the DHI program in the mth month 

the number of herds in the Kth category and jth 
month enrolled on the DHI program 


the number of herds in the kth category and jth 
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month not enrolled on the DHI program 
The amount of this milk that was actually delivered to 
plants could be calculated by multiplying this total 
predicted production for each category by an appropriate 
monthly delivery coefficent. Summing over all categories for 


a given month allowed total milk deliveries to be predicted. 


D. The Simplified Models 

To develop these models the breakdown of producers by 
category was observed for June 1979. These 69 producers were 
chosen in accordance with the proportion of herds that each 
category represented for this month. The 69 producers were 
randomly sampled from the 147 producers that were used in 
the development of the complex models. The breakdown of 
producers by category is shown in Table 4.3. 
It was theorized that use of these models should guarantee 
ease of use ines only one equation was necessary for 
predicting as opposed to the five needed with the complex 
models. Regression equations were calculated for the same 
time period for both a Stage-of-Lactation and Cows-in-Milk 
model using the same model specification as in the complex 
determinations. Milk production using this model was 


theorized to be as follows: 


Y(jtm) = N*P(j) + G*n*P(j) 
where: 


Y = the total milk production for the kth category in 
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rr 0 ee 


Category Number 
Fluid 25 
GE 75/77 7 
GE 78/79 11 
1<10200. 14 
1>10200 12 


the j+mth month where j = the base month and 
th Road Oka ROA Lee OR aay 

p = the average predicted production of herds on the 
DHI program in the jth month 

G = a production coefficient for dairy herds outside 
the DHI program in the mth month 

N = the number of herds in the jth month enrolled on 
the DHI program 

n = the number of herds in the jth month not enrolled 


on the DHI program 


This predicted production was then multipled by a delivery 
coefficient. For these models this coefficient was a monthly 


weighted composite of the coefficients determined previously 
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In total the complex and simplified model 


specifications resulted in 432 equations being generated. 
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VY. ANALYTICAL RESULTS 


A. Introduction 

There are three groups of criteria used to judge least 
squares regression results: theoretical, statistical and 
predictive. The theoretical criteria were not important to 
this study since a structural analysis of the coefficients 
was not made. The statistical criteria were used to check 
the hypothesis of a linear model specification. 

The preceding chapter outlined the analytical procedure 
used in the development of an appropriate milk predicting 
model. The choice of the variables used in each model was 
predicated on ‘a priori’ reasoning and the need to use the 
two DHI files as a base. Least squares analysis was used to 
determine the parameters of each models base. Initially the 
hypothesis that the number of cows in the herd and total 
herd production should have a linear relationship was 
tested. Four different model specifications were tested on 
the Complex Stage-of-Lactation model using the Traditional 
Fluid and small Industrial producer(<10200 MSQ) category. 
These two categories represented the two seasonality 
extremes associated with milk production. Regression 
analysis was performed with November’s stages of lactation 
predicting monthly production for each month up to June. 
Model specifications tested were a linear, double 
logarithmic, and two semi-log formulations. One semi-log 


formulation used logged independent variables while the 
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other was log-linear where only the dependent variable was 
transformed. To facilitate these transformations any 
variable observations in the study that were 0 were set 
equal to 1. The following criteria were then used in 
choosing the appropriate model: 

1. the standard error of the estimate24 

2. The coefficient of multiple determination(R2) 

3. Simplicity in use 

Other things being equal, the linear model would be chosen 
because of its relative ease of use and ‘a priori’ 


reasoning. 


B. Results on the Choice of Model 

The results on the choice of model are presented in 
Table 5.1 and Table 5.2. Since the F-tests and t-tests were 
not substantially different among the models tested only 
criteria (1) and (4) above are reflected in the results. The 
linear model was chosen over the double log specification 
mainly on the strength of a lower standard error of 
estimate, its simplicity and ‘a priori’ reasoning. 

Multiple regression analysis was then performed on the 
data unique to each DHI file. Complex Stage-of-Lactation and 
Cows-in-Milk models were developed with monthly indemendent 
variables for each herd regressed with monthly production 
for the months of November 1978 to June 1979. 


24 a small standard error would suggest a better fit of the 
regression line with the data and more confidence could then 
be placed in the results obtained. 
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FLUID MILK CATEGORY 


Month Yield Linear Log-linear Linear-log Dble. log 
Nov. 37874 6156 8588 10601 6380 
(94.8 (82.3) (84.8) (90.4) 
Dec. 39755 6974 9054 1523 Oren 
(94.0) CBDi0) Aes (90.4) 
Jan. 40739 7659 902 4 42550 7159 
(93.1) SBT, (82.3) (89.5) 
Feb. 37483 6914 11484 1753S 1098s 
(93.3) (81.8) Co oy (88.1) 
March 41973 PAY 10470 $3414 8882 
(93.3) (79.5) (80.8) (B55) 
April 41001 7246 10734 12420 9248 
(93.7) (78.0) Bree 7) (83.5) 
May 43250 6664 10223 12356 9007 
(94.6) (S08) (83.0) (84.7) 
June A222 7 6488 9346 11679 8388 
(94.5) (822-0) (83.5) (85.5) 


Month 
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Dec. 


Jan. 


Feb. 


March 


April 


May 


June 


Yield 
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SMALL INDUSTRIAL MILK CATEGORY 


Linear 
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The Stage-of-Lactation model expressed monthly 
production as a function of the number of cows in 6 
different stages. A list of the equations generated for each 
category for these eight months is presented in Appendix 5. 
The Cows-in-Milk model expressed the same basic 
relationships although the five different stages of 
lactation were grouped into one variable. Cows dry was 
retained as an explanatory variable. A list of the equations 
generated for each category for this model is given in 
Appendix 6. 

Two further models were developed on the basis of 
Stage-of-Lactation and Cows-in-Milk. These two models, 
referred to as the simplified models, were developed with 
the idea of determing if milk could be predicted 
irrespective of category. The equations for these two 
developed models are listed in Appendix 7 and 8. 

Although the Cows-in-Milk model was only a modified 
Stage-of-Lactation model, the differences between the two 
with respect to the significance of individual variables and 
the multiple coefficient of determination suggest one might 
have inherent advantages. One such advantage was theorized 
to be that the Stage-of-Lactation model provided more 
information and therefore offered a stronger "signal" of 


future milk production. 
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C. Correlation Analysis 

The presence of multicollinearity problems was checked 
for the four models even though the efficiency of forecasts 
is unaffected by the correlation between the independent 
variables. Forecasts depend only upon the total amount of 
information about the dependent variable contributed by the 
independent variables and is not dependent upon their 
effects taken separately2>5. In this study, any problems 
associated with multicollinearity were anticipated to be 
with the Stage- of-Lactation models since these models had 
the greatest number of independent variables. It was 
theorized that a high degree of association should exist 
among certain variables as they move through time. This is 
due to the fact that cows in one stage at time "t" may be in 
a succeeding stage at time "t+i". For example, cows in the 
41-80 day stage all will be in the 81-150 day stage in two 
months time. To substantiate this, the correlations among 
these variables were observed. Table 5.3 shows the degree of 
association among the November, January and April variables 
in the Simplified Stage-of-Lactation model. The results 
obtained substantiated this hypothesis. For example, cows in 
the 0-40 day stage in November have an increasing degree of 
association with succeeding stages in January until the 
151-305 day stage is reached. A levelling off of this 
association at this stage is expected since November cows 


25David A. Aaker Multivariate Analysis in Marketing: Theory 
and Applications Belmont Calif. Wadsworth Publishing Co. 
1971 
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will not have reached this stage in 60 days time. The 
151-305 day stage in April reveals a high correlation with 
November's 0-40 day stage. The November cows-dry stage also 
appears to have a high degree of association with January’s 
0-40 day stage and less over time with January’s other 
stages. On the other hand, the degree of association of 
November’ s cows-dry stage with April’s stages indicates an 
increase until the 81-150 day stage is reached and a 
decreasing relationship thereon. Both of these relationships 
substantiate the initial hypothesis that any month's dry 
cows could be viewed as a "Slack" variable pointing to 
future milk production. Any November cows that were dry and 
freshened in December or January would be reflected in 


January’s first stage and April’s third stage. 


D. Coefficient Tests of Significance Analysis 

To test the theory that the ‘cows dry’ variable could 
be viewed as a ‘slack’ variable in predicting future months 
milk deliveries the coefficients and their level of 
significance was observed. As expected the level of 
significance for this variable increased through time. This 
can be interpreted as the variable becoming more important 
when predicting milk delivery several months into the 
future. All the models exhibited this significance. In all 
cases the variable associated with the number of 
cows-in-milk was significant at the .01 level. For the 


Stage-of-Lactation models this significance was pronounced 
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with the first four variables. Cows that remained in milk 
after completing the normal 305 day lactation generally were 
not significant at the .10 level. In this respect the latter 
models supplied information that was not possible with the 


Cows-in-Milk models. 


E. Other Results 

The original development of the Stage-of-Lactation and 
Cows-in-MilkK models was done on DHI data. To properly 
ascertain the amount of "bias" or difference between these 
two groups, milk deliveries for the two groups were 
compared. The time period considered for this part of the 
study was July 1978 to June 1979. This table allowed the 
developed models to be used to predict milk production for 
herds on DHI and then using this predicted amount, 
extrapolate over herds in the same category but not on the 
DHI program. The coefficients derived from this study are 
shown in Table 5.5. 

Aggregating this production over categories allowed 
total production in Alberta to be predicted. Since final 
results were expressed as milk delivered to individual 
plants, these production figures had to be adjusted 
downward. 

A random sample of eight producers from each category 
was made and milk production and milk deliveries for the 
same 12 month period was compared. Results are expressed in 


Table 5.6. The average predicted production per herd in each 
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NON DHI PRODUCTION AS A PERCENTAGE OF DHI PRODUCTION 
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DELIVERIES OF MILK AS A PERCENTAGE OF PRODUCTION 
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category for the Complex Stage-of-Lactation and Cows-in- 
Milk models is shown in Table 5.7. 

The total number of herds for each category and the 
proportion of these herds that was also enrolled on the DHI 
program was determined. The difference between total herds 
in a category and the number of herds on DHI for the same 
category allowed calculation of the number of herds within 
category and not on the DHI program. Using herd 
identification numbers from the DHI program and cross 
referencing these numbers with their Milk Board registration 
number allowed separation of herds into the appropriate 
categories. Using the Summary of Milk Producers information 
for May 1978, May 1979 and June 1979, the number of herds 
for each category was estimated. To determine the number of 
producers both on and not on the DHI program for use in the 
simplified models, a summation of herds for each of the 5 


categories was made. Table 5.8 summarizes these findings. 


F. Predicting Milk Deliveries using the Complex Models 

The actual amount of milk delivered for each month over 
the years 1974 to 1979 (Figure 5.1) was obtained from the 
Alberta Dairy Control Board monthly milk and cream 
production data. Using the average predicted production per 
head for each category, milk deliveries were predicted over 
the same time period, November 1978 to June 1979, used in 
the development of the models. Prediction time periods 


ranged from the present month to as many as 7 months into 
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AVERAGE MONTHLY PREDICTED PRODUCTION BY CATEGORY 
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the future. Table 5.10 indicates the percent error 
associated with each predicted period. With a few exceptions 
predicted production was low. Estimated milk deliveries 
displayed the greatest deviation from actual deliveries when 
present month forecasts were made but tended to be more 
reliable when future deliveries were predicted. In large 
part estimates tended to vary the most when November and 
December time periods were used as the base for prediction. 
This was due to a large number of herds who were in the 
Industrial milk category these two months moving into the 
GE78/79 category in January. This movement influenced future 
predicted deliveries since these herds were now 
representative of a category with a higher average 
production per herd. The predicted deliveries over the 7 
months were determined to have an average 2.0% error. On a 
per cow basis this amounted to an error in estimation of 
2160 cows-in-milk if each cow averaged 4091 kgs. per year. 
An example of the calculations involved in predicting 
March’s milk deliveries using November as the base month is 
shown below. The general November equation for predicting 


milk for a category was as follows: 


Y = D*N*P + D*n*G*P 


where: 
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PERCENT ERROR PER TIME PERIOD USING THE COMPLEX MODELS 
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total milk delivered for March 

delivery coefficient for a given category 
March 

number of herds enrolled on the DHI program 
a given category in November 


average predicted production for a herd in a 


given category for March 


n = the 


number of herds not enrolled on the DHI 


program for a given category in November 


G = the 
the 


production coefficient for dairy herds outside 


DHI program in March 


Therefore the predicted deliveries for March for each 


category was as follows: 


Fluid category: 
£3952) 0274) 444973) fey ( $952) inog666).03%1) (47973) 


20,821,709 kgs. 


GE 75/77 category: 
(.944)(74) (30662) + (.944)(.769) (111) (30662) 


42612,635 Keser 


GE 78/79 category: 
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2,834,175 
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I <10200 MSQ category: 
(.953)(38)(11162) + (.953)(.901) (445) (11162) 


= 4,669,227 


I >10200 MSQ category: 
ee 9 (47 15258 es (6839). ( 1.0023) ( 276 ):( 15253) 


= 4,214,754 


Summing over all the categories the total predicted 
production for March was therefore 37,152,500 kgs. The 
actual milk delivered for this month was 37,667,818 kgs. The 
percent error associated with the prediction was therefore 


1.37% low. 


G. Predicting Milk Deliveries using the Simplified Models 
These two models were based on an aggregation of 

certain data components used in the complex models. These 
aggregations are shown in Table 5.9. Appendix 9 details the 
calculations involved in determining the two composite 
production and delivery coefficients. Using the average herd 
production, a composite production coefficient of .865 to 
compensate for herds not on DHI and a composite delivery 


coefficient of .933 to reflect the actual amount of milk 
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delivered, deliveries of milk over the same time period used 
in developing the models was done. The percent error 
associated with each predicted period is shown in Table 
5.11. Predicted production for these models tended to be 
both above and below the actual amount 75% of the time. The 
average error of estimate was 5.5%, a greater variation than 
that obtained in the Complex models. An example of the 
calculations involved in predicting March's deliveries using 
November as the base month was as follows: 


Y=) D#N*Pe+ D*n*#G*P 


where: 

Y = the total milk delivered for March 

D = the composite delivery coefficient of all herds 
in March 

N = the number of herds enrolled on the DHI program 
in November 

p = the average predicted production for a herd 
in March 

n = the number of herds not enrolled on the DHI 
program in November 

G = the composite production coefficient for a herd 


outside the DHI program in March 
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Therefore we have: 


oss) 522) (2506 ies 933) (865) (1249) (25087) 


= 40,199,032 kgs. 


The actual amount of milk delivered in March was 37,667,818 


kgs. Thus the predicted value was 6.72% high. 
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VI. Summary, Conclusions, and Implications 


A. Summary 

This study focused on two aspects of predicting milk 
deliveries in Alberta. Firstly, it attempted to ascertain 
the feasibility of developing confident milk predicting 
models using the DHI program as a base and secondly, if so, 
which method of approach is was most appropriate for 
obtaining accurate results. The period of study was November 
1978 to June 1979 and four models were tested. Two of these 
models were grouped and called complex models and were 
fitted to monthly production data of the Dairy Herd 
Improvement program by using ordinary least squares (OLS) 
techniques. A linear formulation was used in developing the 
models. 

Another two models, referred to as simplified models, 
differed from the former models with respect to herd 
producer designation. The complex models consisted of herds 
separated into five producer categories on the basis of past 
production trends and assumptions about future production. 
The simplified models consisted of a random sampling of 
herds irrespective of category. Each of these model 
groupings were differentiated into two types of models on 
the basis of the data used in the study. One specification 
was developed solely from monthly cow data with cows 
separated into 5 stages of milk as well as cows dry. The 


other specification used monthly herd data with the number 
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of cows-in-milk and number of cows-dry as the only 
independent variables. 

To properly account for those herds outside of the 
Dairy Herd Improvement program and in the same producer 
category, deliveries between a randomly sampled number of 
herds both on and not on the program was made. A further 
comparison of the amount of milk delivered relative to the 
amount produced was made for each category from a random 


sampling of eight producers on the DHI program. 


B. Conclusions 


"Its easy to predict the future. 


The hard part is being right." 26 


The data examined and analyzed, while subject to 
certain limitations of accuracy, consistency, and 
reliability does enable some conclusions to be drawn. 

1. It is possible to predict milk production and milk 
deliveries in Alberta using least squares technique and 
data from the Dairy Herd Improvement program as a base. 
The so-designated complex models were able to predict 
production over the same period of time used in the 
development of the models with an average error of only 


2%. Using the simplified model versions predicted 


26 Harry Browne, New Profits from the Monetary Crisis New 
York: St Martens Press 1979 p42 
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deliveries had an error of 5.5%. Both model 
specifications were able to predict in the face of an 
overall increasing trend in milk deliveries and a 
significant delivery falloff in the month of February. 
The significance that the cows-dry coefficient achieves 
in predicting milk deliveries into the future 
substantiates the hypothesis that the number of cows-dry 
in the base month can be viewed as a slack variable. 
That is, these cows eventually freshen and take on an 
importance in future milk production. 

Considering the above conclusions, it is suggested that 
the Complex Cows-in-MilK model be used to develop the 
study further. This model simplifies the 
Stage-of-Lactation model into just two variables, 
cows-in-milk and cows-dry. Data needed for testing of 
the Stage-of-Lactation models is not currently available 
and a further 18 month period would be needed to fully 
develop and confidently test these models. The 
Stage-of-Lactation model does provide a much stronger 
signal of a farm manager’s herd intentions since milking 
cows enter more specific stages than just a broad stage 
of cows-in-milk. For DHI inter-herd movement this is 
important since milking cows that enter the herd from 
the purchase of another DHI herd would enter the herd by 
being placed in those stages of the lactation curve that 
correspond to the number of days in milk for each new 


cow. If a preponderence of these cows were to fall into 
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one stage their effect might more readily explain future 
deliveries. The Cows-in-Milk model does not take this 
information into account. 

On the other hand herd enlargement through herd 
purchases probably takes place among herds not enrolled 
on the DHI program. Accurate days-in-milk for each cow 
is not available and the simple designation of whether 
each cow is dry or in-milk is the maximum information 
available. This is particularly true when small 
industrial producers and cream producers are involved 
since so few of them participate on the DHI program. The 
dynamic nature of this segment of the dairy industry is 
due to many such producers upgrading or expanding their 
operations to meet graduated entrant requirements, 
expanding to become larger industrial milk shippers or 
being sold to other dairy producers who are expanding. 
Because this latter movement would be the most 
prevalent, the Complex Cows-in-Milk model is preferred. 
The simple information of cows-in-milk and cows-dry 
meets the needs of the Cows-in-Milk model unlike the 
Stage-of-Lactation model. In addition, the Complex 
version with its producer categories is preferred to the 
Simplified version since it produced more accurate 
results. 

Increased accuracy and confidence in the model is 
obtained if producers are initially partitioned into 


specific producer categories. The main advantage is that 
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the standard error of the estimate is much less when 
producers are initially categorized. 

5. A breakdown of the herds presently enrolled on the DHI 
program suggests that in the last year new enrolments on 
the program have mainly come from the cream category. It 
appears that those milk producers interested in the 


program are already participants of it. 


C. Limitations of the Study 

Although predictive accuracy was better than first 
expected considering the large variability among producing 
herds, certain limitations to the study preclude definitive 
conclusions being drawn. Firstly, the degree of accuracy 
possessed by the models was compared over the same Hine 
period used in the development of the models. This should 
result in more accurate results than those expected if the 
models were tested outside of this time period. Secondly, 
the models were developed and compared over the relatively 
short time period of eight months. This was made necessary 
since data for development of the Stage-of-Lactation models 
could only be collected on an “as is" monthly basis and 
would have required 18 months lead time if models 
encompassing a six month future time period were desired. 

Thirdly, the volatility of the industry as it rapidly 
changes in response to a more open graduated entrant policy 
makes it difficult to segment herds into appropriate 


producer categories. This aspect could be the major reason 
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for the models tending to predict low. The simplified models 
overcome this problem to some extent since only the total 
number of herds producing in the base month is needed rather 
than the requirement that the movement of herds among cream, 
industrial and GE 78/79 categories be clearly monitored. A 
monthly determination of the number of producers by category 
would be necessary for model reliability. 

Fourthly, the amount of milk delivered as a proportion 
of herd production (Table 5.7) raised questions of the 
validity of some of the data. In many instances, herd 
production for a selected month was less than the milk 
delivered by that herd for the same month. If this had 
occurred on an occasional basis, it could be explained as 
herd expansion through purchases not being reflected in the 
DHI production data but being felt in the amount of 
deliveries. However the large number of occurrences suggest 
either a fundamental error is present in the manner in which 
some milk figures are determined or producers are not 
reporting all the cows in their herd. Since deliveries form 
the basis of payment to producers and deliveries occur 
approximately every 3 days for a producer, the delivery 
figures would appear to be accurate. 

Lastly, although the induetnial sector is producing 
less of the total production in the province, the absence of 
many DHI herds in the industrial sector placed limitations 
on the sampling involved in developing the models. This 


resulted in an inability to compare DHI deliveries with non 
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DHI herds for the same producer category. 


D. Directions For Future Research 

Two major areas provide possibilities for future study. 
Firstly, now that milk deliveries have been predicted based 
on estimating average production per herd, accurately 
determining the number of dairy cows in the Province would 
prove useful. Currently these numbers are obtained from 
Statistics Canada studies and are updated every five years. 
This work could prove beneficial if the number of cows in 
milk and the number of dairy heifers in the province also 
formed a part of the study. In both these areas the work 
itself would supply valuable information as to the direction 
of the industry (expansion or contraction). 

A second area of possible study would be an attempt to 
explain the difference in monthly milk deliveries by the 
larger industrial milk producers (Table 5.6). Deliveries as 
a percentage of production for these producers differ 
substantially from those of the other producers. Barring 
sample error, the differences here may reflect management 
considerations not characteristic of the other producer 
categories. An in depth study of sampled producers from each 
of the five categories would test this hypothesis. 

Insofar as this particular study is concerned it would 
serve a useful purpose if it was completed so that monthly 
usage was possible. In main, this completion is enhanced due 


to the simplicity with which the four milk predictive models 
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were developed. It was assumed that milk production could be 
predicted by simply using the number of cows-in-milk as the 
only important variable. In a larger supply response study, 
and over a longer time period, this variable would be 
coupled with other explanatory variables chosen on an ‘a 
priori’ basis. Considering the Complex models predict easily 
over the short term and appear to possess a high potential 
for accuracy, it is doubtful if the benefits associated with 
adding variables to improve the accuracy would not be more 
than offset by the increased costs involved in collecting 
further formation on a monthly basis so that the model could 
be used regularly. 

The ease with which this program can form an integral 
part of the information system of Alberta Agriculture can be 
readily seen. Computer capability exists as a result of the 
DHI program and the monthly routine of obtaining reliable 
future delivery totals is quite simple once the Siren 
herds, equations and pertinent tables form part of the 
system. In the interests of ease of use base month 
predictions should probably be limited to four, eight, and 
twelve month t ime periods rather than the monthly ones 
involved in this study. To allow further development of 
these models a listing of all the herds that participated in 
the study is given in Appendix 10. 

In large part the development of milk predictive models 
such as outlined in this study allows for the accurate 


assessment of the dairy industries response to factors that 
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influence decisions at the farm level. 
that no procedure is in place to monitor these changes, 
use of a milk predicting model has validity and the 


completion of this work has much merit. 
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VIII. APPENDIX 1 


UNIT RETURNS-SETTING MECHANISM FOR 


Canadian Dairy Commission 


i ___, 


Ey 


Direct Subsidy to Milk 


Producers 


INDUSTRIAL MILK 


Indexed via Returns 
Adjustment Formula to: 
- Consumer Price Index 35% 
- Cash Input Price Index 45 % 
- Judgment Factors 20 % 


Canadian Dairy Commission 
Gross Target Return on 
Industrial Milk 


Canadian Dairy Commission 
Offer-to-Purchase Scheme 


Support Price of 
Skim Milk Powder 


8.2320 kg Skim Milk 
Powder per hectolitre 
of Milk 


| Support Price 
of Butter 


4.3218 kg Butter per 
hectolitre of Milk 


Market Price Guarantee to 
Processors in Milk Equivalent | 
Terms (3.6 kg Fat per hL) 


Less 


Estimated Producer 
Market Return 


Assumed Processing 
Piant Margin 
EF 


Gross Target Return on 
Industrial Milk 
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Tx. APPENDIX 2 
Alberta Dairy Development Program Expanded2? 

The recently announced change to the Alberta Dairy 
Development Program is intended to provide cream producers 
and industrial milk producers with an opportunity to change 
to fluid milk production and to provide fluid milk producers 
with an opportunity to expand their dairy enterprise. It is 
also intended to help anyone interested in having a dairy 
enterprise, to get started. 

Assistance under this program is given in the form of 
interest rebates on loans obtained from lending agencies 
that have a guarantee agreement with the Alberta 
Agricultural Development Corporation and which come under 
the corporation’s Farm Development Loan Program. The amount 
of the interest rebate is $70 per $1000 borrowed for a 
period of three years. 

The change in the program involves the maximum loan 
that will qualify for the interest rebate. In the case of an 
individual, it is now $30,000 instead of $12,000. In the 
case of a partnership or corporation it is now $50,000 
instead of $18,000 for two people and $24,000 for more than 
two people. Under the present program, a partnership or 
corporation must consist of more than one farmer or farm 
family to qualify for an interest rebate. 

An applicant who received an interest rebate on a loan 
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of $12,000 in the past would be eligible for an interest 

rebate on the difference between $12,000 and $30,000. 

Similarly, a partnership or corporation that received an 

interest rebate in the past on a loan of $18,000 or $24,000 

would be eligible for an interest rebate on the difference 
between either of these figures and $50,000. 

The rebate can be claimed after a complete annual 
installment of principal and interest has been paid by the 
borrower. If however, the borrower is in debt to the 
corporation either directly or through guaranteed advance 
loans being in arrears, it may apply proceeds from the 
rebate against the debt. 

Following are examples of the types of projects that 
would qualify for an interest rebate. 

1. The construction, alteration, reparation or extension of 
a building for dairying purposes on a farm already owned 
or being purchased. 

2. Permanent improvements required by the dairy enterprise. 

3. The purchase of dairy equipment includes equipment used 
in final processing, handling and storage of feed. 
Examples include: forage wagons for "OQ" grazing, silos, 
forage conveying systems in and to dairy production 
facilities, mixer mills and grain rollers. 

4. Items directly relating to manure handling systems 
including: manure spreaders, pumps and even small 
tractors with blades or loaders to be used exclusively 


in manure handling. 
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Hay storage sheds and facilities for raising replacement 
cattle and caring for dry dairy cows. 

The purchase of milk quota. 

Other purposes related to a dairy operation that have 
the approval of the Alberta Agricultural Development 


Corporation. 
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KX. APPENDIX 3 


INSTRUCTIONS TO THE DAIRY HEROD IMPROVEMENT 
PROGRAM USER 


To ensure a meaningful herd report, the user has an obligation to 
provide all data as required to maintain complete records for his herd. 


THE REPORTS ISSUEO BY THE DHJI. PROGRAM CAN ONLY 
REFLECT THE ACCURACY OF THE ORIGINAL RECORDINGS. Computer 
programs will check the incoming entries to see if they are acceptable. If 
not, the remarks column on the herd report will signal the unacceptable 
recording. Incorrect entries are not directly altered during computer 
processing. The herd owner must submit the corrections to maintain 
error-free records. 


HOW TO MAKE CHANGES 


If any changes are necessary they may be made on the BARN SHEET 
by striking out the incorrect information and writing in the correction 
immediately below. An explanation of the change may be mede in the 
“REMARKS” column. 


IMPORTANT NOTICE! 


Since the cow inventory (index) numbers are the computer's only 
means of identifying the individual cows in your herd, their numbers must 
nat be changed without prior consultation with, and the approval of, the 
Dairy Herd Improvement Data Control Co-ordinator. 


EXPLANATION OF HERD REPORT 


. HERO 1D. - A unique identification used to identify the herd to the 


system. The first two digits indicate the year of enrollment of 
the herd on the D.H.). Program. 


REGION & DISTRICT - A code used to identify the D.A. district in which 
the herd is located. 


OWNER - Name and address to which the herd report will be mailed. 
TEST DATE - Day on which the sampling and recording was completed 
MAILING DATE - Date the computer processed the herd. 


TESTING PLAN - Owner-sampler or Supervised. Information provided on 
both programs is identical. The supervised testing plan 
provides an authenticated record of production. 


BATCH - For Data Control use only. 
REG. - Indicates if 50%, or more, of the herd is registered. 
COW IDENTIFICATION 


BREED - refer to “BREED CODES” on reverse side of this page 
A registered cow is indicated by a ‘1° preceding the breed code. 


REGISTRATION OR NIP 
OR EARTAG NUMBER - AI registered animals of each breed are 
¥ permanently identified by this number as shown 
On their registration certificates. Grades must be 
identified by NIP or eartag numbers. 


BIRTHDATE - Although the cow’s birthdate is kept on record, only the 
month and year of birth are printed. 


LAST CALVING DATE - The date on which the cow began her current 
lactation. 


LAST REPORTED BREEDING DATE-The most recent breeding date 
reported for the cow. 


TEST DAY INFORMATION 


. STATUS - refer to “ENTRY STATUS CODES” on reverse side of this page. 


MILK - The total of the P.M. and A.M. milk weights reported on the barn 
sheet on the test day. 


_ % FAT -The actual butterfat test recorded on the barn sheet for the test 


day. 


- GRAIN FED - The amount of grain mixture or concentrates fed to each 


cow.on test day as recorded on the barn sheet. 
COW MANAGEMENT MEMOS 
RECOMMENDED GRAIN DAILY 


The daily grain recommendation is based on the voluntary intake of 
forage and the qualities of the forage and the grain fed The foilowing 
factors are taken into consideration in the calculation of the daily grain 
requirements on an individual cow basis: 

- bodyweight 

+ age 

- stage of pregnancy 

- milk production on last test day 

+ average fat test during current lactation 

- voluntary intake of forage 

= quality of grain fed 

The daily grain recommendation is the calculated amount of grain 
required to meet the cow's individual needs based on information reported 
on the barn sheet. It can only be as accurate as the reporting ard may 
have to be adjusted by the dairyman after careful observation of his cows 
The minimum recommendation is 2 kg. 
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COW MANAGEMENT MEMOS 


Status Change 
Code. Required 

DB due to breed A predicted date, 60 days beyond the last 
calving date is printed for all cows in milk 
between 36 and 74 days 

OPEN cows not bred The word OPEN is printed for all cows in 
"milk more than 74 days and not reported 
bred. 

EG pregnancy check All cows in calf between 27 and 70 days 
will have a pregnancy check date printed 
which is 42 days beyond the breeding 
date 

eye) due to dry For all cows milking and in caif between 


170 ard 210 days, a due to dry date will 
be printed that is 222 days beyond the 
breeding date. 


LF lead feeding For all cows incalf between 211 and 229 
- days, a lead feeding date, 264 days 

beyond the breeding date will be printed 

oc due to calf 278 days beyond the breeding date for all 
cows in calf between 230 and 275 days 

(Average gestation period for Holsteins). 


COW INVENTORY NUMBER - The number assigned to the cow for 
computer identification. 


BARN-NAME OR NUMBER - The name or number used by the owner to 
identify each cow. Up to 7 letters or 
numbers (including spaces) are allowed 


LACTATION NUMBER~ The actual lactation number (or number of 
calvings) of each cow, regardless of whether the 
animal was previously on test or not. 


BODY WEIGHT (TENTHS) - The weight of each cow expressed in units of 
10 ts estimated from the heart girth measure- 
ment the first month on test and re-estimated 
the first test after each calving. 


CURRENT LACTATION 


Current lactation indicates the number of days a cow has been IN 
MILK, DRY, or IN CALF since the most recent calving date, as well as the 
amount of milk and fat produced since that date. 


For cows that are considered “Record Commence” the production 
totals will include only the production accumulated from the date they are 
entered on test. Predicted production cannot be made for these animals 

a 


BREED CLASS AVERAGES 


B.C.A.'s are indexes of the cow’s production expressed as percentages 
of the 305 day Breed Class Standard, for all cows of the breed calving at 
the same age. For a cow that has been in milk for more than 45 days and 
less than 305 days, estimated 8. C.A.’s based on current lactation totals 
are given. This estimate assumes that the cow will be in milk fora full 305 
day lactation. The deviation from herd average is the difference between 
the cow’s B.C.A. milk and fat index and the Rolling Herd Average B.C.A. 
milk and fat index of all cows completing lactations in the herd within the 
fast 12 months 

Breed Class Averages are an evaluation of the production ability of a 
cow adjusted for age and month of calving B.C.A.’s enable the dairyman 
to compare each cow's performance agarnst the performance of ai! cows 
im all herds. 

The deviation from the Herd Average allows the dairyman to make a 
comparison of a cow’s performance with that of other cows milking in the 
same herd, under the same conditions. 

BCA. values provide a reliable means. through this deviation from 
the herd average, to be used in culling low producers from the herd Herd 
replacements should be kept from superior cows. 


PRODUCTION AVERAGES 


Production averages are expressed for two different time periods: 
- Herd average for the most recent test 
- Herd average for the past twelve months 


These figures are beneficial to the herd owner to assess where his 
herd ranks in relation to other herds. Production Averages, % Days in Milk 
and Calving Interval can be used to detect problems in genetic ability, in 
feeding and in other herd management practices. 

The daily average production of milk per cow this month is calculated 
on the basis of the number of cow days in milk. 


~ AVERAGE NUMBER OF COWS is the average number of cows in the 
herd (both milking 2nd dry) for the report period after adjustment for part- 
time cows entering and leaving the herd. It is computed by dividing the 
number of cow days on test this month by the days in the report period 


% DAYS IN MILK is the number of cows in milk during the report 
pericd expressed as a percentage of the average number of cows (previous 
column). Number of cows in milk ts adjusted for part-tume cows calving, 
going dry or leaving the herd during the period 
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INSTRUCTIONS TO ASSIST YOU IN COMPLETING 
THE BARN SHEET E 


“CURRENT TEST INFORMATION 


This area should be completed for each cow in milk. 


SAMPLE No. - The bottle containing the milk sample taken from each caw 
is identified by this number. 


MILK - Each milk weight for the day is recorded to the nearest tenth. Milk 
weights less than 10 receive a zero as the first digit. e.g. 8.5 = 
08.5 


FAT % - Will be recorded in the Centralized Milk Testing Laboratory. 
GRAIN FEO - Record the total amount of concentrate fed per day. 


STATUS CODES - If a change in status occurred for any cow since last 
day of test, an entry indicating the day of occurrence 
and the corresponding code to the event is expected 
Space for two changes in the same period is provided, 
eg. 6 16 04=dry April 16th, and 1 03 O5 = calved 
May 3rd. + 

For further explanation of status changes refer to 
“STATUS CHANGE CODES”. 


CURRENT TEST DATE - The date the A.M. weights were taken. 


TEST DAY MILK YIELD - Weight of milk shipped to dairy plant on test day 
(bulk tank reading for two milkings) 


NEW COW INFORMATION 


Each cow to be entered on the program must be recorded with the 
information required. All data must be complete before a permanent 
record can be created. 

Sire registration number is a must for all new cows as it plays an 
important role in the sire evaluation aspects of the Dairy Herd 
Improvement Program. 

Oam indentification will insure the automatic updating of pedigree 
information. Special forms which were provided on a semi-annual basis 
for updating of cow and dam pedigree information should no longer be 
required. 

For entry status, please refer to “ENTRY STATUS CODES” 


For new cows requiring a2 status change, please refer to "STATUS 
CHANGE CODE”. 


PRICE AND FEED DATA (based on last test data) 


FORAGE - The quality and amount of forage fed is keyed to voluntary 
intake. ¥. 


MILK PRICING - Net price of milk. average fat test and fat differential 
received by the dairyman. 
This information provides a current guide to cost of production and 
profit through: 
- total feed cost 
-- income on milk produced 
- income over feed costs. 


. FORAGE CODES (voluntary intake) 


1. Fair 12% of body weight 

2. Good 2% of body weight 

3. Excellent 2'%% or more of body weight 

4. Energy intake estimated from laboratory results. 


GRAIN CODES (value based on energy level} 
1. Excellent 3. Fair 
2. Good 4. Poor 


SPECIAL MANAGEMENT MESSAGES 


With every herd report two lists of Special Management messages 
8re produced: 

- herd management messages 

- Cow management messages 

The purpose of the Special Management messages is three-fold 

1. to alert Data Control staff of potential errors related to production, 
breeding and feeding 

2. to slert the dairymen to individual cow or herd situations which 
fequire adjustments in management practices. 

3. to identify herd or cow production levels that are exceptionally high 
end compliment excellent management practices. 

The Special Management messages are intended to focus attention 
on information that is already reported on the herd reports; they are useful 
to the D_H_1 Technic:an in assisting herd owners in the interpretation of 
HERO REPORTS. 
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APPENDIX 5 


THE COMPLEX STAGE OF LACTATION MODEL 


The following pages outline the coefficients to be used 


in predicting milk with the Complex Stage-of-Lactation 


model. The functional relationship hypothesized is as 


follows: 
Y = the 
X1 = the 
X2 = the 
X3 = the 
X4 = the 
X5 = the 
X6 = the 


predicted 


number 
number 


number 


‘number 


number 


number 


OF 
of 
of 
of 
of 


of 


production for a given herd 


cows-in-milk 
cows-in-milk 
cows-in-milk 
cows-in-milk 


cowsS-in-milkK 


from 0 - 41 days 
from 41 - 80 days 
from .6 7) 150 days 
from, 151) >) 305 days 
for 306+ days 


dry cows in the herd 


By multiplying predicted production by the number of herds 


on DHI and adding this to the predicted production for herds 


not enrolled on the DHI program total milk production for a 


given category can be determined. Multiplying this figure by 


a weighted delivery coefficient calculates total milk 


deliveries for the category. Summing over categories allows 


total milk deliveries to be calculated. 
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Kuli APPENDIX 6 
THE COMPLEX COWS IN MILK MODEL 


The following pages outline the equations to be used in 
predicting milk with the Complex Cows-in-Milk model. Data 
for this model’s development was obtained from the DHI Cow 
Master file. The functional relationship in the equation is 
as follows: 

Y = the predicted production for a given herd 
X1 = the number of cows-in-milk for a given herd 


X2 = the number of dry cows in the herd 


The herd categories are as follows: 


FEUID® =" fradi tional 7flufeetproducers 
GE75/77"%= Graduated entrants from 1975°to 1977 
GE78/79 - Graduated entrants from 1978 to 1979 
1<i0200 - Industrial producers (less than 10200 MSQ) 
I>10200 - Industrial producers (more than 10200 MSQ) 


By summing over individual herds for a given category total 
milk production for the category can be determined. 


Aggregating categories gives total monthly production. 
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TABLE A6.1 


PRODUCER 
CATEGORY 


FLUID 


GETS / 77 


GE78/793 


1<10200 


I1>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING NOVEMBER MILK WITH NOVEMBER‘ S COWS 


ESTIMATED COEFFICIENTS! 


INTERCEPT al x2 R2 DW2 
-4776. 1 64.2644 114.9% .93 2.464 
(38.3) (143.0) 

-11700.9 875.8*** 281.3  .37 1.79++ 
P9568) (34219) 

- 1996.6 G27 Seam 174.8 89 . 2.424 
(68.2) (172.9) 

-2582.0 A75.2Qe** 273.2% .61 2.14+ 
(113.0) (340.6) 

4586.3 296. 5%** = - 59 1.76+ 
(59.8) 


HHH HH HEH HE HE HE PE POP EHP PHP HHP EPH HF HP HE HPP PPE PEPE PEE 


2K AE OK 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.2 


PRODUCER 
CATEGORY 


FLUID 


GE/5/77 


GE78/79 


1<10200 


1>10200 


Se 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING DECEMBER MILK WITH NOVEMBER’ S COWS 


ESTIMATED COEFFICIENTS! 
INTERCEPT yal x2 R2 W2 


“SOW 4 <1 Se Set ce I 5 
(43.4) (162. 


OF oo arg ae 
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We Oy (C180: 


“Devore AOA or te UEOee Ts 254 tees 
Kiso) (41073) 
4815.1 I) iat id baie 8) wow Weoote 


(78.4) (200.4) 
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standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE AS.3 


PRODUCER 
CATEGORY 


PEO ED 


GELS/ 77 


GE78/79 


I<10200 


I>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING JANUARY MILK WITH NOVEMBER‘ S COWS 


Eo hiMaleED COEFFICIENTS! 


INTERCEPT rel 12 R2 ss DW2 

-5696.3 686.0*** 369.2%* .90 2. 18++ 
(4% .%) (178.3) 

-8949.9 905.1*** 140.6 .48 1.66++ 
pads. 5), (262028) 

- 1204.0 532.3** 247.6 .82 1.67++ 
(Foes 4902.2) 

-7284.7 531.4*** 993.8%* .56 1.93+ 
(138.4)  (49%.2) 

5091.1 79.044 299.4 67 2.0444 


CH2r 6) (4 85'..7) 
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standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for ser-ial correlation is inconclusive 
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TABLE A6.4 


PRODUCER 
CATEGORY 


FLUID 


GEVSAT 


GE78/79 


1<10200 


I>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING FEBRUARY MILK WITH NOVEMBER’ S COWS 


Eo EMORED CORBPPICTIENTS! 


INTERCEPT Mi ie. R2 — Dw2 

-4474 3 633.5*** 273.2* .90 2.19++ 
R445)  (166:..2%) 

-6079.4 805.0*** 88.0 .54 1.61+ 
“Ute 25) (€19%.8)) 

393.4 444,9*** 316.7* .80 1.74++ 
(At. 7) (187.8) 

-6224.9 467.0*** 944.9*x* .54 2.11+ 
(128.0) (385.9) 

4383.6 260.8%** 232.4 .65 2.20++ 


fooo2y CBT. tt) 


ft HF HEHE HE EH HE HE HP AP HPF HP HPP ER PH HPP TPP HP ETP HPP POPP PPP EPP 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE AG.S 
THE COMPLEX COWS IN MILK MODEL 


PREDICTING MARCH MILK WITH NOVEMBER’ S COWS 


PRODUCER ESTIMATED COEFFICIENTS’ 
CATEGORY INTERCEPT X41 oO) Ro Dye 
FLUID -4729.3 736.3*** 167.6 .90 2.20++ 
(49.7) (185.6) 
GE75/77 -6480.2 878.8*** 202.8 .60 1.56+ 
(15%.7) (195.5) 
GE78/79 2540.0 A64.4*** 345.7 .76 2.19++ 
(83.9) (212.9) 
1<10200 -5983.7 510.4%** 966.3% .52 2.15+ 
(144.2) (434.6) 
1>10200 5117.9 258.6%** 357.9% .62 1.91++ 
(86.5). (224.2) 


HHHHHHHHHAAAHRAR RAR EARHAH RRR RRARA RAR RRR RRA RARE 


: - standard errors given in parentheses. 


- Durbin Watson statistic at the 0.01 level 
KK - coefficient significant at the 0.01 level 
* - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
++ - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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TABLE A6.6 


PRODUCER 
CATEGORY 


FLUID 


GE75/77 


GE78/79 


1<10200 


I>10200 


THE COMPLEX COWS IN MILK MODEL 


INTERCEPT 


-4631.6 


725 ome 


3209 «2 


Bee aiaks 


3970. 0 


4) 


ai) 


ai) 


0) 
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ie 


vhteE 166. 
(165. 


ase te clio. 
(146. 


Giese evel 
(247. 


~okee 7126. 
(434. 


comet as 200k, 
C20: 


PREDICTING APRIL MILK WITH NOVEMBER‘ S COWS 


PomIMareD COEPEICIENTS' 


Ré DWe 
.91 2. Ofias 
69 1.44+ 
sO) cok 
47 1.93+ 
ne) 1.66% 


HHBHHHHRHARERRA ARERR EH EE EE 


2K 2K OK 


- standard errors given in parentheses. 


Durbin Watson statistic at the 


coefficient significant at the 


coefficient significant at the 


coefficient significant at the 


G20 
On Ont 
0.05 level 


level 


leve ] 


0.10 level 


no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.7 


PRODUCER 
CATEGORY 


FLUID 


GE75/77 


GE78/79 


1<10200 


I>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING MAY MILK WITH NOVEMBER’ S COWS 


ESRIMATED COBFFICTENTS' 


INTERCEPT De x2 R2 DW2 
-4203.7 700.8*** 378.2** .94 1,84++ 
(39.3) (146.7) 
-2925.2 735.0%** 403.7***.72 1.514 
(108%) (134.8) 
3086. 1 513.6*** 427.6 .63 2.28++ 


Hi2o ik U32R4) 


-4186.1 BOO .4eee SO0.9% .48 t45% 
CTPisth (SES. 7) 


Ce tiss Aie.8*** 438.8% .70 t. 864+ 
FINS .6) (303.2) 


HHH HHA EE ER Ho PE EO PE 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.8 


PRODUCER 
CATEGORY 


FLUID 


GE15/77 


GE78/79 


1<10200 


1>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING JUNE MILK WITH NOVEMBER'S COWS 


INTERCEPT 


- 3455.4 


“396674 


3228.4 


-4865.8 


263533 
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674. 
(36: 


036% 
(141. 


514. 
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a ae 
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918: 
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Chedsus Cley 4 
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Bret 13435 
Rasy 


9) 


5 ) 


UE s 3 0a7 
(3265 


Drees 
(526): 


Axxe 387. 
(300. 


7) 


va 


6) 


EOUIVALED CORFEIGIENTS 


R2 DW 2 

Bx*x*x 94 1.69++ 
9) 

Ox .62 te e2Ts 
5) 

g 163 252677 
7) 

Ox 51 1.94+ 
4) 

2 EES: 2207 7* 
e) 


HH Ht EE FP EEO PP PHP PPE PEPE EPP 


3K OK K 


_ 


standard errors given in parentheses. 


Durbin Watson statistic at 


coefficient significant at 


coefficient significant at 


coefficient significant at 


the 0.01 level 


the 0.01 level 
the 0.05 level 


the 0.10 level 


no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.9 


PRODUCER 
CATEGORY 


FLUID 
GEIST 7 
GE78/79 
T<10200 


T>10200 
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THE COMPLEX COWS IN MILK MODEL 


PREDICTING DECEMBER MILK WITH DECEMBER’ S COWS 


ESI IMATED COEFFICIENTS! 


INTERCEPT X14 we Ree Dw. 

-5493.6 737 .8%e* 964.0" .93 2.28++ 
(39.9) (148.4) 

-11741.5 1000. 2**« - Jay AL 5+ 
[226 .2) 

- 1568.4 555.2*** 171.1  .83 2.03++ 
(re Oy a Oi Bay ates 

-3020.7 469.0*** 420.7* .50 1.65+ 
(188.7), (303.2) 

3747.6 377.6*** -74.2 .74 2.06++ 
oon ATT a 


HHH AEH EH eH EP oP FP a Pt PEE PEE PEEP 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.10 


PRODUCER 
CATEGORY 


FLUID 


GE75/77 


GE78/79 


1<10200 


I>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING JANUARY MILK WITH DECEMBER’ S COWS 


ESTIMATED COEFFICIENTS! 


INTERCEPT et x2 R2 

-5562.9 711.Q*** 255.2* (91 
(46.7) (173.8) 

-8391.7 956.8%** -56.0 .53 
(193.5) (280.6) 

-1107.6 524.6*** 308.1* .81 
(64.0) (165). 1) 

-4278.9 480.4*** 620.6% .51 
C14iGr, St. CSSit., 4) 

4416.0 Said. eee 170, 74 


6 
Ve on WLAO.2) 


DW? 
2 Nk 


1.36 


1.84++ 


bot 


2 vont 


HHPHHHHAHEHRARRA REPRE AEA RR ARERR A ERA RAARARA REE 


KK OK 


standard errors given in parentheses. 


Durbin Watson statistic at the 0.01 


coefficient significant at the 0.01 


level 


level 


coefficient significant at the 0.05 level 


coefficient significant at the 0.10 level 


no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE AG. 11 


PRODUCER 
CATEGORY 


FLUID 


GE75/77 


GE78/79 


1<10200 


I>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING FEBRUARY MILK WITH DECEMBER’ S COWS 


ESTIMATED COEFFICIENTS! 


INTERCEPT al Le R2— DW2 

-4411.8 626.0*** 310.5* .90 2.17++4 
(44.3) (165.0) 

-4966.5 B22.i*et -48,9 .56 1.294 
(155.7) (225.8) 

143.6 430.9"** 412.3 .78 1.94+ 
(61.6) (158.4) 

- 3538.0 406.0*** 639.1* .49 2.05+ 
(135.2) (305.9) 

4052.5 286.9%** 184.8 .70 2.28+ 


P03) (487.7) 


HHHHHHHHHHREARAR ARERR HARE AR ARERR REAR RE 


2K OK AE 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.12 


PRODUCER 
CATEGORY 


FLUID 


GEIS/ 77 


GE78/79 


1<10200 


I>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING MARCH MILK WITH DECEMBER’S COWS 


ESEIMATED COEFRICIENTS! 


INTERCEPT ma X2 R2 _Dw2 

-4775.4 711.3*** 290.4% 190 2,.24++4 
(49.6) (184.6) 

-5655.7 886.5*** 105.7 .60 1.31+ 
56%) 1227.8) 

2384.1 A55.6*** A14.1** .74 2.3144 
2.7) ) (9871.8) 

-3513.9 465.7*** 646.2* .50 2.10+ 
(148.2) (335.2) 

4726.4 280.4*** 327.2* .68 1.91++ 


(90.7) (220.4) 


HHAHHHRREAAEAREAEARAARA PERERA ERR RARR REAR RARER ERRRARRE 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE AG. 13 


THE COMPLEX COWS IN MILK MODEL 
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PREDICTING APRIL MILK WITH DECEMBER’ S COWS 


PRODUCER 
CATEGORY 
FLUID 


GE75/77 


GE78/79 


1<10200 


I>10200 


ESTIMATED COEFFICIENTS" 


ie) 


Roe re 20). 
(1670.4 


saree 282), 
(214, 


4) 


5 ) 


OTE ATS), 
(318. 


Sree 252 . 


ae 


INTERCEPT Val 
-4672.1 695 
(43, 

-3850.9 750 
(115 

3294.5 486 
(8S. 

-1804.5 463 
(140 

5438.5 270 
(84 


a) 


Le 


s0Tee 260 
tie2. 


(204 


4) 


ne 
92 


.68 


.68 


250 


.66 


Aas at 


PaG2r 


TaG3t 


HHHRHRHAR RHA R RAHA RRA RHR AAA AKA A HA HRR ARAB 


1 a 


Durbin Watson statistic at 
coefficient significant at 
coefficient significant at 


coefficient significant at 


the 0.01 
the 0.01 


standard errors given in parentheses. 


leve 1 


level 


the 0.05 level 


the 0.10 level 


no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.14 


PRODUCER 
CATEGORY 


FLUID 


GETS /77 


GE78/79 


1<10200 


1>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING MAY MILK WITH DECEMBER’ S COWS 


ESTIMATED COEFFICIENTS 


INTERCEPT a ie R22 Dw2 
-4165.9 688.3***« 440.5***,.94 1,84++ 
(3%.9) (140.9) 
-3585.2 754.3%%* 367.7** .73 1.184 
(105.5) (153.0) 
3434.0 570.9%** 228.1 .61 2.31++ 
(709.5) (264.8) 
-3028.3 583.4*** 548.6% .53 1,.59+ 
(162.6) (367.8) 
671.4 538.5%** 317.8 .77 1.82++ 
(118.7) (288.4) 


ee 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TREBLE AG. 15 
THE COMPLEX COWS IN MILK MODEL 


PREDICTING JUNE MILK WITH DECEMBER’ S COWS 


PRODUCER ESTIMATED COEFFICIENTS! 
CATEGORY INTERCEPT X4 Ne R2 DW2 
FLUID - 3368.5 679.4%** 346.3%* .95 1.64++ 

(S5y 5) 132 43) 

GE75/77 -4367.6 829 .5*** 254.4 .65 1.54+ 
(186.5) (197.9) 

GE78/79 3578.4 553. 4*** 255.9 61 2.30++ 
(10628) (2806...1) 

1<10200 -3754.2 647. 9*#* 513°9+ 57 1. 80+ 
(163,68) > °( 370.5) 


I>10200 2O0CH 2 S50. Eee 356.9 8 2,044 
(ines) 4029082) 


HEHEHE HEHRRAARARAAERRRARARAARARARERRRRRAR RARER RRR AE 


: - standard errors given in parentheses. 


2 - Durbin Watson statistic at the 0.01 level 
> > ok - coefficient significant at the 0.01 level 
* - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
++ - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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AAOLE AG. 16 
THE COMPLEX COWS IN MILK MODEL 


PREDICTING JANUARY MILK WITH JANUARY’ S COWS 


PRODUCER ESTIMATED COEFFICIENTS? 
CATEGORY INTERCEPT Ma x2 R2 DW2 
FLUID -5017.1 7oa.2%%4-187.0 .94 .2.08++ 

(43.8) (972.2) 

GE75/77 -4563.50  639.8*** 432.6%**.89  1.38+ 
(100.4) (46.3) 

GE78/79 - 1997.1 538.5*** 357.8%**,.82 1,96++ 
(56.0) (134.3) 

1<10200 -4490'3° 463,7*** 395.2) .62 4.97++ 
(115.1) (336.6) 

1>10200 4122.4 407.0***-114.2 .81 2.06++ 


(62.9) (152.2) 
HHEHHEHRAAARHRARRRERARRRAARAEAAAAERRRARARRAEARAR ERE 


: - standard errors given in parentheses. 


cs - Durbin Watson statistic at the 0.01 level 
7K - coefficient significant at the 0.01 level 
> - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
++ - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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TABLE AG. 17 
THE COMPLEX COWS IN MILK MODEL 


PREDICTING FEBRUARY MILK WITH JANUARY’ S COWS 


PRODUCER ESTIMATED COEFFICIENTS! 
CATEGORY INTERCEPT X14 x2 R2 DW2 
FLUID -3997.7 693.8x** -48.6 92 2.31++ 

CaO Atti. 3} 

GE75/77 -1919.8 574.2*** 336. 3***. 86 1.55+ 
(93.2) (43.0) 

GE78/79 -894.4 451.2e** 437. 3***. 80 1.97++ 
Cbe.6)-  ¢126.5) 

1<10200 -988.7 373.24*% 497.4% .57 2.05++ 


Cia Oy) S26. 2) 


I>10200 3838.0 Std, OF#* 106. 4 ve Ae et 
(74.4) (179.9) 


HHHHHHEREHAREHRARARRARAR ARR ARRARAARARREA RARE ARH ER RE 


: - standard errors given in parentheses. 


2 - Durbin Watson statistic at the 0.01 level 
7 > 2k - coefficient significant at the 0.01 level 
2k = coetricient significant at the 0.05 level 
7 - coefficient significant at the 0.10 level 
++ - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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TABLE A6.18 


PRODUCER 
CATEGORY 
FLUID 


GETS (TT 


GE78/79 


1<10200 


T>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING MARCH MILK WITH JANUARY’ S COWS 


ESFIMATED COERFICIENTS' 


INTERCEPT na xO R2 = Dw2 

-4335.0 762.Q%** = - 292) 240+ 
(29.4) 

-426.1 632.7*** 351.0***,.89 1.42+ 
(90.6) (41.8) 

865.2 473. 7*** 492. 4xx*.78 2.3144 
(60.9) (145.6) 

-573.2 426.4"** 474.6 .57 2.154 
(122.2) (358.7%) 

4429.1 315.3*** 234.0 .69 1.83++ 
(00.0) (000.0) 


ee 


1 — 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.19 
THE COMPLEX COWS IN MILK MODEL 


PREDICTING APRIL MILK WITH JANUARY’ S COWS 


PRODUCER ESTIMATED COEFFICIENTS? 
CATEGORY INTERCEPT X41 X2 R2 Dw2 
FLUID -4947.5 735.9*** 43,7 93 2,25++ 
(450) (174,27) 
GE75/77 1526.30 604,.5*** 249, 0*** , 87 1.32+ 
(79.3) (36.6) 
GE78/79 1655.0 492 .5%**x* 409.9** .72 2,33++ 
(7023) (1686.4) 
1<10200 633.7 ALA Saxe O77 5+ .6if 1.96+ 
(Ino). (325.8) 
1>10200 5075.6 328.0*** 74,3 .69 1.65+ 


(80.8) (195.5) 
HHHAHRERRAAAEAAAAAHARARBRRRA RRR R ARBRE AAARR RRR RRARRRRE 


1 - standard errors given in parentheses. 


2 - Durbin Watson statistic at the 0.01 level 
OK - coefficient significant at the 0.01 level 
2K - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
++ - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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TABLE A6.20 


PRODUCER 
CATEGORY 


FLUID 


GE/5/77 


GE78/79 


1<10200 


1>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING MAY MILK WITH JANUARY’ S COWS 


ESLIMATED COEFEICIENTS! 


INTERCEPT al X2 R2 Dw2 
-3815.9 730.0*** 209.1% .95 1.91++ 
(35.8) (149.5) 
2399.7 644.0%** 215. 1%"*,.84 1.55+ 
(89.4) (41.3) 
1748.6 563.4*** 408.8% .64 2.33++ 
(95.1) (228.0) 
159.5 558.2*** 274.6 .64 1.83+ 
KIZPAOh) (372.4) 
483.4 573.1*** 222.3 .78 1.83++ 
(tse) (2602) 


HHHRHEHREARERERARARRRRARAARARAREAREARRRRREARRR ARREARS 


1 


standard errors given in parentheses. 


Durbin Watson statistic at the 0.01 level 


coefficient significant at the 0.01 level 


coefficient significant-at the 0.05 level 


coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.21 


PRODUCER 
CATEGORY 


FLUID 


GE75/77 


GE78/79 


1<10200 


T>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING JUNE MILK WITH JANUARY’ S COWS 


EOP iMateD COEFFICIENTS! 


INTERCEPT ak X2 R22 DW2 
-2944.3 723.8%** 93.3  .96 1.58+ 
(82.9) (137.5) 

2128.2 612.9x** 318.3***.90 1.48+ 
(79.2) (36.6) 

1965.7 548.2*** 421.4% .63 2.34++ 
(94.9) (227.3) 

-107.1 596.5*** 245.7 .65 2.00++ 
(133.3) (389.7) 

2188.4 528.2*** 499.7* .78 2.08++ 


(1te.4) (288.9) 


HHHHHHAH HAHAHA AERRR RARER RRA RRR RARER RRR RAR RRR ARR 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 


152 


S205 WstM WT meee 1319095 


2) Yaad MRE LM ahi eit 


no 
bt 
Lr (Mato 12 498GD 433G) Cala AM! 2d 
yd a oe %, ix 


a 
te SE) NRRED 


2 eats aye «nee [One ae SSS 3 
(@ =i 


ar (SOs 
(1D + end B68 1. eae7 
fe See). ey ; 
5 UNS wana) See ' /Ot- 
ebe) . (f;8et) ie 
és , . ay ’ Cop iris AOA cgi 
(@,.88$) Petty 


5049 VERE NNSA ANN 


JSS ey ah Tayi 21S ve 7 i 
level 104th te ob aitate noo ew nidnud - - 
wei $0.0. sett is ingeryrnghe nai Shao = *~ 

vel 20.0 se ing heumegd thrygh a” rotor tees - 
Laval ‘O)70ane $3 lta tigy 2° nein eag. + 
ira st: Noh etSi10 lay neg heel 


eds 


avigul agin? . at anand dtd Tat tow sot feat odd = a pL + 


TABLE A6.22 


PRODUCER 
CATEGORY 


FLUID 


GE75/77 


GE78/79 


1<10200 


I>10200 


iss 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING FEBRUARY MILK WITH FEBRUARY’ S COWS 


Eo LMaED) COEFF IICTENTS"' 


INTERCEPT Xai 
-4044.8 702 
(31 

-5128.8 642. 
(85. 

1965.2 468. 
(60, 

-1705.9 537. 
(48. 

-10.2 601 
(136%. 


OD 


ore 162}. 
nO) Qi 4. 


Gres 330 


7) (38. 


Tee 130). 
tei (TAR. 


Axxe 190. 
One hosp 


Oren 43 Oo}. 
9) (432. 


R2 DW2 
5 94 2 1a ak 
4) 
o*#* BO 1. 36* 
3) 
4 . 86 24, 20H 
i) 
ome . 80 2 Oe 
3) 
1 Oe 1.624 
3) 


HAHAHAHA HR EE PH EE EE 


standard errors given in parentheses. 


Durbin Watson statistic at the 0.01 level 


coefficient significant at the 0.01 level 


coefficient significant at the 0.05 level 


coefficient significant at the 0.10 level 


no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE "AOO23 


PRODUCER 
CATEGORY 


FLUID 


CEO 11 


GE78/79 


1<10200 


I>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING MARCH MILK WITH FEBRUARY’ S COWS 


EOUiMAlED COEFEICTENTS 4 


INTERCEPT x4 eo R2 = DW2 

-4440.7 758.0*** -48.3  .93 2.20++ 
(38.4) (165.0) 

-2848.0 675.5*** 353.4e**.89 1,294 
85-9) (38.3) 

2364.1 A94.7*** -80.9 .80 1.83++ 
(80.9) (154.5) 

436.9 506.6*** 373.0** .82 2.73+ 
(64.9) (134.2) 

281.9 668.2***-492.1 .68 1.59+ 


(148.9) (470.2) 


HHEHRHHAHRHEHRHRARRRARARRRERRARARARARARAR RARE 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.24 
THE COMPLEX COWS IN MILK MODEL 


PREDICTING APRIL MILK WITH FEBRUARY’ S COWS 


PRODUCER PouIMAneD COBPEICIENTS' 
CATEGORY INTERCEPT Kt X2 Re DW2 
FLUID -4381.5 720.4*+* (55.4 .94 220044 

(34.5) (148.3) 

GEID/77 822.0 094 .6*** 265.5*** . 86 teags 
(86.3) LOR. 7) 

GE78/79 30H SiS 400). OF **- 137-2 8 coo ht 
(20.2) (147.4) 

1<10200 1200.5 A60.2*** 498.6** .76 2.47+ 
9.3) (964.2) 

1>10200 1048.9 650.5+**-506.4 aii 130% 


(134.3) (424.0) 
HHAHAHHHAHHAAARRRAH RAR RAR RARER RRR RARER RRR 


U - standard errors given in parentheses. 


2 - Durbin Watson statistic at the 0.01 level 
2 ak 2 - coefficient significant at the 0.01 level 
* - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
++ - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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TABLE A6. 25 
THE COMPLEX COWS IN MILK MODEL 


PREDICTING MAY MILK WITH FEBRUARY’ S COWS 


PRODUCER ESTIMATED COEFFICIENTS! 
CATEGORY INTERCEPT X4 Xo R2_Dw2 
FLUID -3837.2 727. 5*** 156.0 .96 1.70++ 

(Sn te C125. 2) 

GE75/77 2982.4 595. 2xx* 243,2e**.80 1.63+ 
(t0623). (44.7) 

GE78/79 - 1549.7 727. 5eee = 81  2.03++ 
(84.3) 

1<10200 1061.5 504. 1*** 597.7** .67 2.37+ 
(096% ~(226).8) 

1>10200 386. 1 756.0***-457.6 .71 1.714 


(Pi@2e.Ger GOs; 9) 
i 


: -- standard errors given in parentheses. 


2 - Durbin Watson statistic at the 0.01 level 
* 2K - coefficient significant at the 0.01 level 
* k - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
++ - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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TABLE AG.26 


PRODUCER 
CATEGORY 
FLUID 


GE/5/77 


GE78/79 


1<10200 


I>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING JUNE MILK WITH FEBRUARY’ S COWS 


EooMoLeD COREELCIENTS! 


INTERCEPT iG. X2 R22 Dw2 
-2949.7 Pita tert 250,39" 964. Tage 
(2605) (1328) 

3034.1 560. 1*** 346.9***.88  1.58+ 
(87.0) (38.8) 

Das? 673.0*** 189.6 .80 2.10++ 
(128.4) (245.0) 

1005.7 545.3*** 460.0* .68  2.49+ 
(10747) (2230) 

-87.8 802.7***-484,9 .72 1.73+ 
(169.0) (533.4) 


HHRHRAHRRARARRRR ARERR AR ARRRR RRA RRR RAR RARE RRR 


1 = 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.27 


PRODUCER 
CATEGORY 


PLU TD 


GE75/77 


GE78/79 


I<10200 


1>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING MARCH MILK WITH MARCH’ S COWS 


EO LIMATED CORFETCLENTS.! 


INTERCEPT Mal Ke R2 Dw2 
-4708.9 826. 4xx*-362. dex 94 2,39++ 
(95.5) A164 
-445 3 611.0*** 369, 8%**,.92 1,394 
(S47) (70.2) 
3370.9 513.6***-258.0* .87 1.78++ 
(60.2) (149.0) 
-512.8 552.2*** 310.7** .86 2.89+ 
(6519) (108.7) 
-645.2 poenorre A005 | 1A. 4 es 
(146.5) (361.4) 


HHHHHHHH RAHA RAR ARR RR RRR RRR RRA Raa HR RRR ARR RRA AA 


standard errors given in parentheses. 


Durbin Watson statistic at the 0.01 level 


the 0.01 
05 


10 


coefficient significant at leve | 


coefficient significant at the 0. level 


coefficient significant at the 0. leve] 


no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE AO. 28 
THE COMPLEX COWS IN MILK MODEL 


PREDICTING APRIL MILK WITH MARCH’ S COWS 


PRODUCER ESTIMATED COEFFICIENTS’ 
CATEGORY INTERCEPT “al X2 R2 Dw2 
FLUID -4623.9 770.3***-168.6 .94 2.21++ 

(36.2) (154.4) 

GE75/77 2402.6 555. 6*** 267.0%**.87 1.444 
(nay 2) (36.8) 

GE78/79 4536.4 396.0*** -60.8 .76 1.48+ 
(72.8) (180.9) 

1<10200 575.3 AQ8.4*** 431.2***.80 2.624 
(79,8) (129.8) 

1>10200 -26.0 619.7***-128.0 .68 1.63+ 


(149.3) (368.4) 
HHEHAHHEHHRAAREAARARE ARERR REAR RAR RRA R AAR 


: - standard errors given in parentheses. 


NO 
| 


Durbin Watson statistic at the 0.01 level 


2 2K - coefficient significant at the 0.01 level 
> - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
a - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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TABLE A6.29 


PRODUCER 
CATEGORY 


FLUID 


GEIL RT 


GE78/79 


I<10200 


I1>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING MAY MILK WITH MARCH’ S COWS 


ESPIMATED, COEFRICTENTS' 


INTERCEPT i X2 R2 DwW2 
-4036.7 763.3*%** 9 - 96 1.79++ 
(20.8) 
3662.5 579, 3%*** 236. 1***, 2 1.76++ 
(87.6) (43.3) 
-239.0 657.8*** 55.4 82 1.93++ 
(109.1) (269.9) 
288.4 556.7*** 498.7** .72 2,52+ 
(iS) (0184.6) 
-893.1 656.5*** 166.6 .66 1.66+ 


(183.6) (453.0) 


HHHPHHHHHAAAHEARA RARER AR RRR RA ARARAARAREEE  RR  RP 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is -inconclusive 
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TABLE 6.30 


PRODUCER 
CATEGORY 


FLUID 


GE75/77 


GE78/79 


I<10200 


I>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING JUNE MILK WITH MARCH’ S COWS 


bol IMeVED COEFFICIENTS | 


INTERCEPT x) KO R2 — Dw2 
-3228. 1 739.4%** -28.2 (96 1.54+ 
(30.4) (129.6) 

4420.9 527. 9*** 350.8%**.90 1.9444 
(25,59! 487.3) 
1058.7 703.4***  - 84 2.2444 
(89.6) 
-574.9 635.5*** 315.4% .73  2.58++ 
(106.6) (178.2) 
-1440.5 725.6*** 91.8 .68 1.87++ 


(188.2) (464.5) 


HHHHHHHHHHAEAEA HERE RAPP A RRR ARARR ARAB RARER REE 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE AG... 37 


PRODUCER 
CATEGORY 


FLUID 


GES TT 


GE78/79 


I<10200 


I>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING APRIL MILK WITH APRIL'S COWS 


ESRIMATED COEFFICIENTS! 


INTERCEPT X41 EE R2 Dw2 

-5488.9 763.4*** -58.6 .96 1.73++ 
GSO) 4. 3) 

-96.6 600.3*** 270. 3***,.90 1.37+ 
(67.5) (32:..3)) 

-22.0 554. 8*** 287.4** .83 2.67+ 
Gop.6)° (i3se.3} 

303.7 533.6*** 77.7 .69 1.46+ 
CRRA Ok: CSMa Sh 

A79T.3 409.7***-110.9 .83 1.54+ 
(52.6) (128.8) 
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2K AK OK 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.32 


PRODUCER 
CATEGORY 


FLUID 


GE75/77 


GE78/79 


1<10200 


1>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING MAY MILK WITH APRIL’ S COWS 


EolimAleD COEFFICIENTS! 


INTERCEPT X41 12 R2— DW2 
-4607.2 722. 5%** 2AQ.7e* 196 1.57+ 
(S053) (11522) 

399.0 141. 5%** 232. 5%#* 86  1.80++ 
(80.2) (38.4) 
-468.0 614. 7*** 367.2* .76 2.59+ 
(gOsa) (19207) 
-1000.8 586.5*** 402.6 .68 1.47+ 
(136, 4) (382. 5) 
127.8 686.3*** -42.3  .86 1.94++ 


(79.6) (184.4) 
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standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 


VOD 2 we HT EW ALIN YAM + Ung 


~~, 


PAN EIeNG EE ts agrawite? 
eS ovo Sys dk 


f ae) #> ee ae ‘ 
1 peta 


ae aa? . rs +4ng , ww 
& ae 6S, 08) 


pr oc: 23 we Tae weet Sra 
(T.ser) (e204) 


rn) ee 8, 20m prema Joga ae 
(2.280). Cec TRO)” 0) 0 a ee 
as 6 SB wane ady oe a 
(6 bap) 9a 8T) 7 


2 ky SAO ROHR ENON ARAN 
assonirened. HF. ay ana 
Fowl POD a ts or neitene noataw nt 
Oo ett fe reatooore meta rs 

aveal 29 ease, 46 
\ tevel db,Doaat te 


{rise 


| Vere : 


4ytéeuforsant et reapers) 


RABLE AG?.33 


PRODUCER 
CATEGORY 


FLUID 


GET 5/77, 


GE78/79 


T<10200 


I>10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING JUNE MILK WITH APRIL’S COWS 


EOLIMATED COEFFICIENTS’ 


INTERCEPT X14 X2 Re pwe 
- 3687.9 670.9*** 348.7***,.95 1.564 
C0 Bie RR 
2224.4 572.0*** 351.3%"*,91 2.1344 
(7'3)..0) (35): 63} 
-1085.9 BS. tere 175.4 0H 2.524 
(99.0) (184.4) 
-1464.7 635.3*** 373.9  .70 1.69+ 
(142.0) (395.2) 
1573.6 714.3*** -58.9 .87 2.26+ 


(3020) (185.3) 
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standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0-05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.34 
THE COMPLEX COWS IN MILK MODEL 
PREDICTING MAY MILK WITH MAY'S COWS 
PRODUCER ESTIMATED COEFFICIENTS 
CATEGORY INTERCEPT im 1D R2— DW2 
FLUID -4761.0 750.0*** 106.3  .97 1.38+ 

(26-5) * (97.3) 

GE75/77 -249.9 665.7*** 193.6%**.88 1.74++ 
(76N 0), (9.5) 

GE78/79 - 1068.8 702.1*** -42.1  .80 2.49+ 
(95.9) (276.8) 

1<10200 - 1228.6 DOO. y* . O72 69 1.66+ 
(300.5) (1172.8) 

1>10200 -244.9 668.7*** 131.6 .90 2.54+ 
(64.9) (124.0) 
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standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABEE: AG)..35 


PRODUCER 
CATEGORY 


FLUID 


GE 75/77 


GE78/79 


1<10200 


EA10200 


THE COMPLEX COWS IN MILK MODEL 


PREDICTING JUNE MILK WITH MAY’S COWS 


ESTIMATED COEFFICIENTS' 


INTERCEPT Ma 0 Re pw 
-3785.2 687.9*** 250.7** .96 1.434 
(28.3) (104.0) 

1947.5 581.3%** 328. 5%#* 91 2.2344 
CBG.) (375.3) 
-629.4 703.7***-120.9  .79 2.18++ 
(97.8) (282.3) 
-1429.4 572.0* 468.8 .68 1.95+ 


(S24) ( W264. 1) 


1154.6 647 ..3#** 251.3* .85 PES) baa 
G7iGI Oe (i tS2t.3}} 
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standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE A6.36 
THE COMPLEX COWS IN MILK MODEL 
PREDICTING JUNE MILK WITH JUNE’S COWS 
PRODUCER ESTIMATED COEFFICIENTS? 
CATEGORY INTERCEPT Bat he Ree DW2 
FLUID -3541.2 682.7*** 233.2% .95 1.78++ 
Poi Oue wuitye4| 
GE75/77 1235.2 594.0*«** 326.5***,92 2.15++ 
(65.8) (35.7) 
GE78/79 -353.4 669.3%** = - 79 2.16++ 
(69.9) 
1<10200 - 1303.0 508.3*** 788.5* .71 1.48+ 
(76 Oie de aio. 5) 
1>10200 -1175.3 677.8%** 321.2* .82 2.68+ 


(dit Sebihes « (f88% 9) 
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standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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KIN. APPENDIX 7 
TRE SIMPETRPEDSSTAGE OF LACTATION MODEL 


The following pages outline the coefficients to be used 
in predicting milk with the Simplified Stage-of-Lactation 
model. The functional relationship hypothesized is as 


follows: 


Y = the predicted production for a given herd 

X1 = the number of cows-in-milk from 0 - 41 days 

X2 = the number of cows-in-milk from 41 - 80 days 
X3 = the number of cows-in-milk from 81 - 150 days 
X4 = the number of cows-in-milk from 151 - 305 days 
X5 = the number of cows-in-milk for 306+ days 


X6 = the number of dry cows in the herd 


Herd categories are dispensed with in the use of this model. 
By multiplying predicted production by the number of herds 
on DHI and adding this to the predicted production for herds 
not enrolled on the DHI program total milk production in 
Alberta can be determined. Multiplying this figure by a 
weighted delivery coefficient calculates total milk 


deliveries. 
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KV APPENDIX & 
THE SSIVPEEE IED COWS IN MILK MODEL 


The following pages outline the equations to be used in 
predicting milk with the Simplified Cows-in-Milk model. Data 
for this model’s development was obtained from the DHI Cow 
Master file. The functional relationship in the equation is 


as follows: 


Y = the predicted production for a given DHI herd 


X1 = the number of cows-in-milk for a given DHI herd 


X2 = the number of dry cows in the DHI herd 


Herd categories are dispensed with in the use of this model. 
By multiplying this predicted production by the number of 
herds on DHI and adding this to the predicted production for 
herds not enrolled on the DHI program total milk production 
in Alberta can be determined. Multiplying this figure by a 
weighted delivery coefficient calculates total milk 


deliveries. 
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TABLE A8.1 


THE SIMPLIFIED COWS IN MILK MODEL FOR NOVEMBER 


PREDICTED 


MONTH 
NOVEMBER 


DECEMBER 


JANUARY 


FEBRUARY 


ESTIMATED COEFFICIENTS 
INTERCEPT ot Xo Re) pe 


“BIS229 182. 07%* 323.2 ee 220 eF 
(4.8) (275.6) 


“$HO25E6 foe, ot*t* B09, Geen. 60 2 oars 
(93.9) (274.8) 


=I206982 162.6*** 772.5***,84 2,.25++ 
(68.6) (254.5) 


“1034915 66a. 2*** 695.3 eo) ee ete 
PSG. G 7 toe) 


HHHHHEHERRHAHEAR RRA RARE RRARRARER RARER AERA ERR RAR AREF 


1 = 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant at the 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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THE SIMPLIFIED COWS IN MILK MODEL FOR NOVEMBER 


PREDICTED 


MONTH 
MARCH 


APRIL 


MAY 


JUNE 


ESLIMAIED COEFFICIENTS ' 


INTERCEPT Xd 
-10677.6 757. 
(64, 

-9297.7 7%2 
(54, 

-7071.8 706. 
(54 

-5033.1 706. 
(56. 


X2 Ro Owe 
Bee 25-744", 65 2a lr 
GS) (238.6) 
pete OUAs eth Oo 2a het 
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ares COOL irre 8G" Daze + 
2) W909) 
WeCr OU ort oo) 2e207F 
BS} 240.9} 
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1 = 


Durbin Watson statistic at 
coefficient significant at 
coefficient significant at 


coefficient significant at 


standard errors given in parentheses. 


the 0.071 level 


the 0.01 level 
the 0.05 level 


the 0.10 level 


no serial correlation is present 


the test for serial correlation is inconclusive 
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FABLE AG .'S 


THE SIMPLIFIED COWS IN MILK MODEL FOR DECEMBER 


PREDICTED ESTIMATED COEFFICIENTS! 
MONTH INTERCEPT at X2 Re.) we 

DECEMBER -10889.1 B10. Fees 407 3209.81 2.23++ 
(77.4) (278.9) 

UANUARY -12185.7 784. 1*** 721,.6%**,.85 2.20++ 
(id. 7) — (258.9) 

FEBRUARY -10531.6 691.4"** 715.9%. 86 2.20++ 
fee Gon) (22970) 
68.5*** 726.7 .86 2.18++ 


MARCH -10868.8 7 
(6656). 9(240..7) 


Ht Hf Hf HH HE HE HP PE oP HPP HP PPP HP EP PPP PEP HEE EE 


L - standard errors given in parentheses. 


“ - Durbin Watson statistic at the 0.01 level 
7 - coefficient significant at the 0.01 level 
* - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
++ - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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TABLE A8.4 


THEE SEMPE LE PEDSCOWSS INSMILKY MODEL: FOR" DECEMBER 


PREDICTED ESTIMATED COEFFICIENTS! 
MONTH INTERCEPT pal ye R22 DwW2 
APRIL -9430.4 729. 4xxx 660.4***.89 2.18++ 
(55.6) (200.4) 
MAY -7151.8 727.7*** 636.2***.90  2.19++ 
(56.61) --(490.9) 
UUNE -5051.9 720.9*e* 484.1** .86 2.21++ 


(5924) (243.9) 
HHHHHEHRERPAREARRHRARE REAR APA RRR RRR RRA ARR REE AF 


! - standard errors given in parentheses. 


2 - Durbin Watson statistic at the 0.01 level 
2 - coefficient significant at the 0.01 level 
* > - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
++ - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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ABLE ABs. 5 


THE SIMPLIFIED COWS IN MILK MODEL FOR JANUARY 


PREDICTED ESTIMATED COEFFICIENTS! 

MONTH INTERCEPT ra X2 R2 DW2 

UANUARY -9884.27 776. 5*** 457, 4x** 94  1,96++ 
(30.8) (43.2) 

FEBRUARY -8256.3 711.3*** 360.6***.93 1.9444 
(30.4) (42.5) 

MARCH -8448. 1 787.5% 360.6%**.92  1.87++ 
(SAO 4%.) 


HHHHHARRHRRARHRAERRHE RRR REAR RAR ARR RRR A RRR 


! - standard errors given in parentheses. 


4 - Durbin Watson statistic at the 0.01 level 
26 a - coefficient significant at the 0.01 level 
7k - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
+ - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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TABLE A8.6 


THE SIMPLIFIED COWS IN MILK MODEL FOR JANUARY 


PREDICTED ESTIMATED COEFFICIENTS! 

MONTH INTERCEPT KH i R2 »=Dw2 

APRIL -7047.1 768,.Q*** 235 .Qxx*,.92 1.8744 
(31.3) (4878) 

MAY -A769.1 772.7*** 196.0***,.92  1.96++ 
(20.6) | 142.8) 

JUNE -3150.3 703.8*** 315.9%* .92 2.06++ 
(29.8) Aner |) 


tH Hf i HF HF it HF HEH HF HOHE HEHE EHP HE HEHE PPE EOP EE PP 


i - standard errors given in parentheses. 


4 - Durbin Watson statistic at the 0.01 level 
2 2k 2 - coefficient significant at the 0.01 level 
* - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
++ - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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TABLE A8.7 


THE SIMPLIFIED COWS IN MILK MODEL FOR FEBRUARY 


PREDICTED ESTIMATED COEFFICIENTS! 

MONTH INTERCEPT Net X2 R2  -_DW2 

FEBRUARY -7653.8 688. 1*** 337.8%**,.94 2.0244 
(24.6) 43723) 

MARCH -7687.2 756. 7*** 344.7%**.94  1,.98++ 
(25. Gi) (42.2) 

APRIL -6055.9 728. 8%*# 234. 1"**.94 2.01++ 
(26.1) (39.6) 

MAY -3605.1 724,.24e* 207.5%**.93 2.04++ 
(26.8)  (40...2) 

JUNE -2307.0 665. 1*** 325,.3"* .93 2. 12++ 
b257 53 (38.7) 


HHHHHHHHHHH HHH HAHAH HHH A ARH RRA AA AH 


: - standard errors given in parentheses. 


2 - Durbin Watson statistic at the 0.01 level 
2k - coefficient significant at the 0.01 level 
7 > - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
++ - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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THE SIMPLIFIED COWS IN MILK MODEL FOR MARCH 


PREDICTED ESTIMATED COEFFICIENTS’ 
MONTH INTERCEPT Sl a R2 
MARCH -8555.3 791. 1%#x 326, deex 94 
(27.6) (46-0) 

APRIL -6916.1 759, 4exx 225 Qexx 94 
(25.6) © (37.0) 

MAY -4478.9 754, 2x 202. 1** .94 
(26.9), 138!0) 

UUNE -3149.9 697. 5*** 309.3** .94 
024.53) (35 22) 


DW? 
1.0344 


1.974 


eg One 


Qe Gas 


HHA HHH RRR RRR RRA RAR HRA RAR HAHAHA ARR RAR A 


i - standard errors given in parentheses. 


2 - Durbin Watson statistic at the 0.01 
xxx  - coefficient significant at the 0.01 
* - coefficient significant at the 0.05 
* - coefficient significant at the 0.10 
++ - no serial correlation is present 
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level 


+ - the test for serial correlation is inconclusive 
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TABLE A8.9 


THE SIMPLIFIED COWS IN MILK MODEL FOR APRIL 


PREDICTED ESTIMATED COEFFICIENTS! 

MONTH INTERCEPT a Le R22 Dw2 

APRIL -6794. 1 749. 4xxx 244 54**.95 1.9744 
M24 tN > Sr) 

MAY -4288.3 741.24%% 295 3xee.94 2.01++ 
(25,6 (3.7) 

JUNE -2895.1 682.4*** 332.5%* 94 2.1344 
(2427 nl BS650) 


Hee HHR HAR HRA RRR AMHR He ARH AKA AAR HR RRR A 


} - standard errors given in parentheses. 


“ - Durbin Watson statistic at the 0.01 level 
> - coefficient significant at the 0.01 level 
* - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
++ - no serial correlation is present 


+ - the test for serial correlation is inconclusive 
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TABLE AS.10 


PREDICTED 


MONTH 
MAY 


JUNE 


THE SIMPLIFIED COWS IN MILK MODEL FOR MAY 


PolmlManeD OEFPICTENTS! 


INTERCEPT Val x2 R22 DW2 
-4020.0 735.7*** 181.6%**,95 2,08++ 
(93.5) 436-9) 
-2452.5 664. 3%*** 319, 9ee*,94 2.1244 
F250) | e972) 


HHEAHHEHEHRRRARAHRARR HRA ARR PARRA RRR RR RR ARR REAR 


1 = 


2K OK K = 


standard errors given in parentheses. 
Durbin Watson statistic at the 0.01 level 
coefficient significant at the 0.01 level 
coefficient significant ete ne 0.05 level 
coefficient significant at the 0.10 level 
no serial correlation is present 


the test for serial correlation is inconclusive 
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TABLE AS. 11 


PREDICTED ESTIMATED COEFFICIENTS!’ 
THE SIMPLIFIED COWS IN MILK MODEL FOR JUNE 


PREDICTED ESTIMATED COEFFICIENTS! 
MONTH INTERCEPT al me R2 Dw? 
JUNE -3031.9 674.2*** 331.4** .94 2.05++ 
Cane) (40.8) 


it HH Ht HO HF HE HF HF HE HEHE HF EPP HEP HPP APE OE HP PP OPP HHO HOH EPP HPP PP 


, - standard errors given in parentheses. 


: - Durbin Watson statistic at the 0.01 level 
#* + - coefficient significant at the 0.01 level 
* - coefficient significant at the 0.05 level 
* - coefficient significant at the 0.10 level 
a - no serial correlation is present 


t - the test for serial correlation is inconclusive 
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COEFFICIENT "GCALCURAT ICONS FOR THEeSIMPLIFIED MODELS 


XVI. APPENDIX -9 


PRODUCTION COERFICIENA: 


Average production 
Average production 
Average production 
Average production 


Average production 


Using the breakdown 


So) (010); Ost. 


G22 tee Os a OSA atte 0( 929) 0135". 


DELIVERY COEFRICIENT 


Average delivery 
Average delivery 
Average delivery 
Average delivery 


Average delivery 


coefficient 
coefficient 
coefficient 
coefficient 


coefficient 


coefficient 
coefficient 
coefficient 
coefficient 


coefficient 


for 
for 
for 
for 


for 


for Fluid category 

for GE75/77 category 
for GE78/79 category 
for 1<10200 category 
for 1>10200 category 


.865 


Fluid category 

GE75/77 category 
GE78/79 category 
1<10200 category 
I>10200 category 


Using the total herd breakdown for June we obtain: 


640(.952) + 161(.946) 
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of non DHI herds in June we obtain: 


+9310 934) te 346.943) one tel. 


Ow. 
per 
1.034 
FePES) 
.998 


998) 


Rie loe 
.946 
.934 
943 
2030 


830) 


es 


‘Ola. Saf + (RAG peer » RED: Ory + (| gse: or + i es 


aoe | al 7 - : aml A 


@ Wowaaan yx | ia 
wt O19} Jae Boa end -auoao 


V ages AG ao iay laphst i900 | not 
Togates OUASTao ‘oF inetal ¥teo0 not to 
mopeine (ORDIRE was Taha ses rior 

= ynageted OOOGIST eh Seiere tabs —hiee 


ae Lan) hh 
‘on 
retao 9m huh 7 eSver Lig mon da rach a 


ri ie 7 
ae | » ' ’ (Mak ‘ a ty 7 (AEG “ppegs ’ (s03,108 0 


+564, 3 ome 4 


iyuignsaay 4 gil 

Mad, .s rrogedec: Ofigl4 ataF inet 91 emg Gia fate 
fo a a 
7: (iogetes TU \ASS oa) dneroeees covileb ® + 
oiauetes ePrergae qo hiptot tae vrevi feb con 
2. = yiogeten 00807 a1 "4 jataie? ago yrevi fob eqete 
pi wisgetigo USOT AT nee shes) Vie weve Tob oyster 
=f | - F 

. 

nistdo sv srivb \o? mabdee rd bas fatot edt at 


“ay 


KI 
tS) ZUR See orn ihe COMPLEX MODELS 


APPENDIX 10 


Following is a list of the DHI herds used in the 


development of the Complex models. By using these herd 


numbers as a base these models can be made to form an 


integral part of the current DHI program of Alberta 


Agriculture. 


2501015 


60D0381 


63D0052 


65D0448 


70D0111 


72D0543 


73D0626 


75D0825 


76D0942 


Fluid Producers 


5400273 5400393 


61D0297 62D0084 


63D0292 63D0478 


66D0337 66D0485 


LODO FID0058 


72D0546 73D0555 


74D0697 74D0744 


75D0853 75D0858 


76D0970 7701004 


56D0034 


62D0142 


64D0042 


67D0286 


71D0097 


73D0561 


74D0746 


75D0860 


CDOs 
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56D0457 


62D0304 


64D0080 


67D0303 


AIDC35 


73D0570 


75D0760 


76D0904 


77D1027 


57D0188 


63D0029 


6500015 


68D0458 


71D0465 


73D0596 


75007977 


76D0906 


7701034 


99D0182 


63D0037 


65D0156 


69D0431 


72D0512 


73D0611 


75D0806 


76D0929 
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56D009 1 


73D0588 


75D0805 


76D0940 


95D0161 


73D0645 


76D0921 


77D0999 


63D0036 


73D0614 


75D0831 


76D0944 


60D0106 


73D0656 


76D0924 


LAD VO 


E 75/77 Producers 


67D0093 68D0070 69D0004 71D0020 73D0569 


73D0640 73D0642 74D0667 74D0675 74D0725 


75D0855 75D0863 75D0881 76D0910 76D0933 


76D0950 78D1148 


GE 78/79 Producers 


60D0451 67D0277 67D0345 68D0085 70D0322 


74D0676 74D0738 7400747 75D0771 76D0894 


76D0954 76D0976 77D0984 77D0990 77D0996 


7701021 77D1075 78Dh120 
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Industria! Producers <10200 Ibs MSQ 


61D0449 67D0087 73D0649 73D0652 74D0685 75D0772 76D0983 


FILDOS88..7 701003) 2700008 77D1030 7701047 77D1056.77D1062 


77D1069 


Industrial Producers >10200 Ibs MSQ 


58D0444 69D0250 71D0026 72D0545 72D0547 73D0602 73D0623 
73D0630 73D0655 74D0678 74D0680 74D0741 75D0770 75D0787 


75D0862 75D0880 76D0905 76D0968 77D1026 
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Lio! POR-nERDSMmEN THE SIMPLIFIED MODELS 


Following is a list of the DHI herds used in the 


development of the Simplified models. 


These herds were 


obtained by randomly samp ling 68 herds used in the 


development of the Complex models. 
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Using these herd numbers 


as a base these models can be made to form an integral part 


of the current DHI program of Alberta Agriculture. 


94D0393 
63D0037 
67D0286 
72D0547 
7300652 
74D0746 
75D0855 
76D0944 
77D0996 


7701047 


Milk Producers 


56D0034 5700188 58D0444 


63D0052 


67D0303 


73D0569 


74D0675 


75D0760 


76D0904 


76D0968 


JAD 1003 


77D 1056 


GsDO0292 


69D0431 


73D0626 


74D0678 


75D0770 


76D0905 


76D0970 


7701004 


77D1062 


63D0478 


70D0322 


73D0630 


74D0680 


T5D077 1 


76D0906 


76D0976 


77D1009 


7701069 


60D0106 


64D0042 


71D0020 


7300642 


74D0685 


75D0772 


76D0910 


76D0983 


FEDI OAS 


77D1075 


61D0297 


66D0337 


71 D035 1 


73D0645 


7400738 


75D0787 


76D0929 


77D0988 


7701018 


78D1120 


61D0449 


67D0087 


72D0545 


73D0649 


74D0741 


75D0806 


76D0942 


77D0990 


77D1030 
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